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Abstract 

We present in this paper the SO^ library, which calculates an analytic decomposition of the Yukawa interactions invariant under 
SO(2N) in terms of an SU(A) basis. We make use of the oscillator expansion formalism, where the SO{2N) spinor representations 
are expressed in terms of creation and annihilation operators of a Grassmann algebra acting on a vacuum state. These noncommu- 
tative operators and their products are simulated in SO^ through the implementation of doubly-linked-list data structures. These 
data structures were determinant to achieve a higher performance in the simplification of large products of creation and annihilation 
operators. We illustrate the use of our library with complete examples of how to decompose Yukawa terms invariant under SO(2A) 
in terms of SU(A) degrees of freedom for N - 2 and 5. We further demonstrate, with an example for SO(4), that higher dimensional 
field-operator terms can also be processed with our library. Finally, we describe the functions available in SO^ that are made to 
simplify the writing of spinors and their interactions specifically for SO(IO) models. 
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1. Introduction 

The orthogonal groups 0{N) and their generalisations have 
played an important role in the construction of modern physics. 
In particular, the special orthogonal groups SO(V) appear nat¬ 
urally in the context of physical systems invariant under rota¬ 
tions, which in turn implies the conservation of the angular mo¬ 
mentum or the determination of the azimuthal quantum number 
for an atomic orbital. The notion of spin used to describe the 
intrinsic angular momentum of particles is another example of 
the importance of special orthogonal groups. Indeed, the spin 
group Spin(V) is a double cover of the special orthogonal group 
SO(N), i.e., Spin(V) is locally isomorphic to SO(N) (see, e.g.. 
Ref. [1]). 

In particle physics, the use of special orthogonal groups 
SO{N) have been very productive in the construction of Grand 
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Unified Theories (GUTs). The original idea of GUT models 
is to embed the Standard Model (SM) gauge group SU(3)c x 
SU(2)l X U(1)y in a larger simple Lie group, so that the three 
SM gauge couplings unify into a unique coupling. The first 
GUT model was proposed by Georgi and Glashow [2] in 1974 
and it introduced SU(5) as the unifying gauge group. The group 
SU(5) has rank 4 as the SM group and the observed fermions 
are grouped in two unique representations 5 and 10, per gener¬ 
ation. 

The possibility of having a GUT model based on the special 
orthogonal group SO(IO) was first accounted by Georgi [3, 4] 
and Fritzsch and Minkowski [5]. The SO(IO) model brought 
new interesting features over SU(5). Each generation of SM 
fermions are accommodated in a unique 16-spinorial represen¬ 
tation of SO(IO) with an additional place for a singlet Weyl 
field, that can be interpreted later as a right-handed neutrino. 
These sterile neutrino states allow naturally to explain the ob¬ 
served oscillations of neutrinos through the Seesaw mecha¬ 
nism [6-10]; giving an extremely light mass to the active neu¬ 
trinos when the sterile neutrino mass is of order of the unifica¬ 
tion scale. The SO(IO) gauge interactions conserve parity thus 
making parity a continuous symmetry. Due to the fact that the 
rank of SO(IO) is 5, there is an extra diagonal generator with 
quantum number B - L as in the left-right symmetric models 
and it is indeed the minimal left-right symmetric GUT model. 
Finally, GUT models based on SO(N), apart from SO(6), turn 
out to be automatically free of gauge anomalies [11]. 

Since the appearance of first SO(IO) GUT model, many 
models based on SO(IO) have been proposed in the litera¬ 
ture (cf. Refs. [12-22] and references therein). In addition, 
other models were implemented within SO{N) unification with 
a rank greater than 5, e.g., SO(12) [23], SO(14) [24], and 
SO(18) [25-27]. SO(18) turns out to be the minimal special 
orthogonal group that accommodates the three SM fermionic 
generations in a unique spinorial representation 256 by choos¬ 
ing properly the breaking chain down to the SM. There are also 
applications of SO{N) as unifying group in the context of mod¬ 
els with extra-dimensions, e.g., SO(IO) in 5D [28, 29], in orb- 
ifold 5D [30] and 6D [30-36]. The group SO(l 1) was also used 
in the context of Randall-Sundrum warped space [37, 38]. 

The breaking of a GUT SO(N) model down to the SM can 
be achieved by different breaking path, with possibly some in¬ 
termediate mass scales. In order to understand the possible 
SO(A^) breaking paths, it is important to identify its maximal 
subgroup (with the same rank as the higher group), so that one 
can express representations of SO(A^) in terms of representa¬ 
tions of the maximal subgroup and therefore understand the 
necessary Higgs sector. In particular, for the group SO(IO) one 
identifies two important maximal subgroups [13, 16], namely 
SU(5) X U(l) and SO(6) x SO(4), which is equivalent to 
SU(4) X SU(2)l X SU(2)r. The first subgroup can be broken into 
the usual SU(5). Instead, the second subgroup can be broken 
into the Pati-Salam model, SU(3) x SU(2 )l x SU(2)r x 11(1)^-/,, 
in which the B - L symmetry of the SM is gauged. It is worth 
to point out that one can also break SO(IO) to the flipped- 
SU(5) [16], where the SM hypercharge is identified with a lin¬ 
ear combination of the diagonal generator of SU(5) with extra 


U(l) generator of SO(IO). 

The purpose of this paper is to introduce the SO^ library 
implemented in the C-n- programming language. The idea be¬ 
hind the conception of SO^ is the decomposition of Yukawa 
interactions invariant under SO(2N) in terms of SU(A^) de¬ 
grees of freedom. This decomposition is particularly useful 
for GUT models based on SO(2A^) that break to an interme¬ 
diate threshold symmetric under SU(A^), since it allows to re¬ 
late the Yukawa couplings in the intermediate theory with the 
GUT Yukawa couplings from the GUT theory, and thus leading 
to predictions. In general, this decomposition can be fastidious 
and error-prone. Our library is meant to simplify this task. 

The SO^ library relies on the oscillator expansion formal¬ 
ism, where the SO(A^) spinor are written in an SU(A^) basis re¬ 
alised through the introduction of creation and annihilation op¬ 
erators of a Grassmann algebra [39]. These operators and their 
algebra are simulated in SO^ by means of doubly-linked lists 
as the appropriate data structure for these problems. This type 
of data structure has higher performance power, since it opti¬ 
mises the memory usage for long chains of operators and the 
data itself in memory do not need to be adjacent. Although the 
SO^ library was projected with the groups SO{2N) in mind, 
it can be easily adapted to the groups SO(2N H- 1) or even to 
other systems where creation and annihilation operators can be 
defined. 

The paper is organised as follows. In the next section, we dis¬ 
cuss the spinorial representations of SO(2A) in a basis in terms 
of the degrees of freedom of SU(A), through creation and an¬ 
nihilation operators defined in a Grassmann algebra. We then 
apply this method to decompose Yukawa interactions invariant 
under SO{2N) in terms of SU(A) interactions. In Section 3, 
we present the general structure of the SO^ library, giving in 
detail the general functions and specific functions for SO(IO). 
In Section 4, we explain the installation of our library and we 
show how to write simple programs. Then in Section 5, we give 
complete examples for computing Yukawa terms in SO(4) and 
SO(IO) with the SO^ library. Finally, we draw our conclu¬ 
sions in Section 6. 

2. The SO(2N) spinor representation 

We review in this section the oscillator expansion tech¬ 
nique [15, 39,40] that is implemented in the SO^ library. This 
technique has been actively explored for explicit computations 
of SO(IO) Yukawa couplings [41^3]. The main idea of this 
technique is to write the two spinor representations of SO(2A) 
in a basis where the spinor components are expressed explic¬ 
itly in terms of SU(A) fields. This is achieved by construct¬ 
ing a Grassmann algebra of creation and annihilation operators. 
One could have used a completely group theoretical approach 
as done in Ref. [44], but the oscillator expansion technique is 
more field theoretical and seems more intuitive to consider the 
case where the breaking of SO(2N) is done down to SU(A). In 
addition there are other methods in the literature [45-48] that 
can be used for computing the SO(2A) invariant couplings, but 
we shall not consider these methods in this paper. 
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We start by introducing the general properties of any special 
orthogonal group SO{N), which are the simple Lie group of all 
orthogonal N x N matrices O such that 

0^0 = OQT =1, (1) 

with the special condition |0| = 1. This group leaves invariant 
the bilinear 

y ^ xiyi + X2y2 + ■ ■ ■ + XNyN, ( 2 ) 

when the N-dimensional vectors x and y transform as 

w, yp —> = Ofjyyy . (3) 

Making an infinitesimal group transformation, the matrix ele¬ 
ments Ofjy can be expanded as 

C)fJV — Sfiy — CO^yAl^y 0{U} ) , (4) 

where is a real antisymmetric tensor, while are N{N - 
l)/2 independent N x A^-matrix generators of SO(N). In the 
vector representation, the generators are hermitian, M^y = M^y, 
and they can be written as 



implying TrM^y = 0, and they satisfy the Lie algebra of S0(N) 
as 


I^A/^y, — i (SijjjAIyp ^pp^vrj ^vri^^p "t" • (6) 

Within the Cartan classification, the Lie algebra associated to 
the group S0(2N + 1) is while to S0{2N) is D^, with N 
being identified as the rank of the algebra. We focus now our 
discussion only on even-dimensional special groups S0(2A^). 
Note that the oscillator expansion technique can also be applied 
to the spinor representation of S0(2N -H 1). 

The spinor representations of S0(2N) can be constructed if 
one introduces a set of matrices {Tp), with fi - 1,.. .N, such 
that 

XiTiJi H- V2r2y2 + ■■■ + x^rj^yM = x^y. (7) 

In order to verify Eq. (7), one must necessarily impose that the 
matrices Tp should obey to: 

{rp,ry) = 2bpy, (8) 

which form a Clifford algebra. It is straightforward to see that 
any ordered product of distinct gamma matrices gives rise to a 
complete set of linearly independent matrices. This fact leads 
to the construction of the so-called spinor representation of 
S0(2N). In Appendix A we give a general proof of the ex¬ 
istence of the matrices Tp. In fact, for any even-dimensional 
Clifford algebra there is only one irreducible representation of 
dimension 2^. Instead of writing explicitly the matrices Tp via 
the 2^ X 2^ generalised Dirac matrices formed from the direct 
product of the Pauli matrices, we write them in terms of a set 


of creation (b^.) and annihilation (bj) operators acting on the 
Hilbert space as 

{b,,b]}^6ij, {bi,bj}^0^{blb]}, (9) 

with i — Each pair bi, b^. of operators can be 

constructed directly from linear combinations of pairs of T- 
matrices as 

bj = + r2j), b] = \{-iT2j-^ + T2j), (10) 

with the inverted relation given by 

T2j-l^-i{bj-b]), Y2j^{bj + b]), (11) 

showing a one-to-one correspondence. General formulae for 
the correspondence between the Clifford and the Grassmann al- 
gebrae are found in Appendix B. 

The advantage of this approach is that one does not need to 
write explicitly the operators bi, b\, one needs only to define 
the vacuum state |0). One defines the Eock vacuum as the vector 
|0) = |0,0,0,..., 0) corresponding to N unoccupied states, which 
is defined by 

f?,-|0> = 0,forallf = 1,...,A. (12) 

One-state vector can then be represented as 

|0> = |0,0,...,1,...,0) , (13) 

where the non-zero entry is at position i. We have just de¬ 
rived the building blocks to construct the spinor representa¬ 
tion of SO{2N) in terms of states obtained from the action of 
the creation operators. Moreover, defining the set of operators 
Tij = b\ bj, it is easy to verify that they satisfy the algebra of 
U(A) as 

[t,7j = SkjTii - SnTkj. (14) 

It is then not surprising to observe that the basis of vectors ob¬ 
tained through the action of products of creation operators on 
the Eock vacuum, 

K - ^l, ’ L < *2 < • • • < ip, (15) 

expands any vector I'P) with coefficients being irreducible fully- 
antisymmetric U(A) tensors, [f/''- '’’. This fact allows us to write 
the spinor representations of SO(2N) in terms of irreducible 
SU(A) tensors. 

2.7. Decomposition into the SU(iV) basis 

The general expression for the spinor representation of 
SO{2N) written in terms of the SU(A) fields is given by 

(16) 

The completely antisymmetric tensors have dimension 

An easy way to compute the dimension of all tensors in 
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Eq. (16) is by noting that it can be read from the Mh-row of the 
Tartaglia’s triangle'. For tensors with large number of indices it 
may be convenient to reduce them with help of their conjugate 
tensors using the Levi-Civita invariant tensor of dimension N, 


and therefore Eq. (16) becomes 
I'E) = |0>^^ + b]\0}f + 


(17) 


2\{N-2)\ 


■ h!" h'!" 


|0) iA/, 


(18) 


{N - 1)! '■> 


■bL. 10)-A. 


+ blbl---bl\0)il,. 

The dimension of the vector space in Eq. (18) is 2^ which is 
in agreement with the dimension of the E-matrices. Within the 
SU(A^) basis, given by the vectors of Eq. (15), any spinor I'T) 
corresponds to a column vector 

= ••• . (19) 

In this spinor representation of dimension 2^, the states'T trans¬ 
form under SO(2A^) as 

'P ^ r = U(cj) T' or I'P') = U(io) I'P) , (20) 

where the unitary transformation U(a>) is given by 

t/(w) = exp ^ j . (21) 

The generators of the spinor representation are constructed 
in terms of the E-matrices as 

V = ^ [r^, Ey] , (22) 

with E^y = E^y and Tr S^y = 0, which guaranties the unitarity 
of U(a>) and |t/(w)| = 1, respectively. It is straightforward to 
verify that S^y satisfies the algebra of SO(2A^) given in Eq. (6). 
It turns out that the spinor representation |'P) with dimension 
2^ given in Eq. (18) is in fact reducible. This fact can easily 
be demonstrate by observing that the product of 2N E-matrices, 
Eo, dehned as 

Eo - i'^rir2...r2^, (23) 


anticommutes with all E^ matrices, but it commutes with E^y 
and therefore splits the spinor T* into two nonequivalent irre¬ 
ducible spinors 'Ey and of dimension 2^^^, given by 

'Ey = i(l+ro)T'. (24) 


'This mathematical representation is also known as the Pascal’s triangle. 
The triangle was already known centuries before in China, India and Iran. 


The projectors |(1 + Eq) are in total analogy with the chiral 
projectors known in the Dirac space (see projector properties in 
Appendix A). The chiral states 'Ey are generated by the action 
of an even number of creation operators on the vacuum state |0), 
while the chiral states T*- are generated by the action of an odd 
number of creation operators. Observing Eq. (18) one obtains 




'r 


, 'E_ = 

<A,7 


<A' 



y * > 


and one can distinguish two cases in SO(2A) : when N is an 
odd integer 'Ey and 'E- are self-conjugate, while when N is even 
'Ey and 'E_ are distinct spinor representations. 

2.2. Yukawa interactions 

In this subsection, we sketch the construction of Yukawa in¬ 
teractions in a generic GUT model ruled by SO{2N) expressed 
in terms of SU(2A) tensor helds. This is particularly useful 
when the group SO(2N) is broken to its SU(2A) subgroup at 
some intermediate scale, since below the breaking scale one 
gets new relations among the Yukawa couplings of the SU(A) 
theory. In what follows, we shall assume that the fermionic 
degrees of freedom belong to irreducible SO(2A) spinor rep¬ 
resentations and the Higgs helds transform as complete anti¬ 
symmetric tensors of SO(2A). Since the GUT gauge group 
commutes with space-time symmetries, it is convenient to write 
the fermionic degrees of freedom in terms of left-handed Weyl 
helds. In order to write the most general Yukawa interaction 
invariant under SO{2N), it is useful to express the transposition 
of U(cxj) in terms of its corresponding inverse matrix U\co) as 

BU\oj)B-\ (26) 

through the matrix B, which has the property 

B-' r;^ B = . (2?) 

In Appendix A, it is shown that such operator B always exists 
and it can be written as 

N 

Y]r^^(-ifY](bk-bl). (28) 

fi=odd k=l 

From the above equation, one sees that the operator B anticom¬ 
mutes with Eq when N is odd, while instead it commutes when 
N is even. When simplifying expressions involving the operator 
B, it turns out to be more convenient to write it with contracted 
indices as 

B = ^ ibj - b]){bk -bl)--- (b, - bl ), (29) 

where is the Levi-Civita antisymmetric tensor with N in¬ 
dices. From Eq. (26), one deduces that the combination 'E'tB 
does transform as 

'ETb ^ 'E'''B = 'E^B U\ca) , (30) 
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and one concludes that for any pair of spinors 'f'l and '{'2 the 
bilinear 'I'|BC ' 'i '2 is invariant under SO(2A^). Due to the fact 
that 'f'l and 'P 2 are assumed as fermionic fields, the presence 
of charge conjugation matrix C ensures that this bilinear is also 
invariant under Lorentz transformations. Furthermore, the ma¬ 
trices transform as 

U(ajfr^U(co)^0,yYy, (31) 

where the matrix O is given in Eq. (4). Using this result one 
can write an SO{2N) invariant Yukawa coupling combining the 
fermion in a spinor representation 'F with the Higgs scalar <p^, 
that transforms like a vector according to Eq. (3): 

(f>^. (32) 

The relation given in Eq. (31) can be generalised to any prod¬ 
uct of E-matrices as 


U{u)^ r^Ev ■ ■ ■ Ep t/(w) = Opp'Ovv' ■ ■ ■ Dpp-Ep-Ev- ■ ■ ■ Ep-. (33) 

This general formula allows us to write the most general gauge- 
invariant Yukawa coupling under SO(2N). In the language of 
creation and annihilation operators one writes 

Yah ('PLi B C-1 E[p.Ep,... Ep„.] IT',,) , (34) 

where the indices k,l - denote the two possible irreducible 
spinors; a and b are flavor indices and the elements Yah of the 
Yukawa matrix; <^p,p2...p,„ with m - 1,...N is a scalar tensor 
fully antisymmetric. Due to the fact that is fully an¬ 

tisymmetric, any E-matrix product in the formula in Eq. (34) 
should also be made fully antisymmetric, and therefore one has 

r[piEp2-Ep,j = ^ ^(-l)'*^rpp(,)Ep2,(2) • • •Ep2,„„), (35) 

■ p 

where the sum runs over the permutations and dP takes 0 for 
even number of permutations and 1 for odd number of permu¬ 
tations. We notice that the general formula given in Eq. (34) 
can also be applied to the case where the SO(2N) spinors are 
taken as scalar fields yielding a pure scalar interaction in the 
scalar potential. In this case, the bracket in Eq. (34) should not 
include the charge conjugation matrix C. 

In some cases, the bracket given in Eq. (34) vanishes auto¬ 
matically. This can be well understood by taking into account 
the general properties of the projectors |(1 + Eq) and the fact 
that Eq commutes with an even number of E-matrix product or 
anticommutes with an odd number. Thus, in the case N + mis 
an odd number, it implies 

(T'* J B C-' Ep.Ep2 • • ■ Ep,„ IT',,) = 0, (36) 

while in the case N + mis an even number, it implies 
(T'* J B C-' Ep,Ep2 • • • Ep,„ IT',,) = 0, (37) 

when k 1. 


Concerning the antisymmetric tensor <^pip 2 --,,„ some com¬ 
ments are in order. There are N distinct fully antisymmetric 
tensors with dimension with 1 < m < N^. Moreover, 

the representation with dimension denoted as Apip 2 
is indeed a reducible representation that can be decomposed 
into two irreducible representations [12, 14]. For N odd into 
two irreducible pairs self-conjugate representations of dimen- 


Sion 


where 




(38) 


(7f) ) 

_ 1 

,0,r-/2v. 

“ 2 

\ 






'pnVn N'.^Py-PnVyVN 


' ' ''J/invn !\I\^P[-PNVI-VnJ 

Instead, for N even, one has the following decomposition 


(39) 


0 0,l,2--,«’ (40) 

where 


0,i-,v 

_ 1 

,0,r-/2v. 

“ 2 

\ 


JfllVi ■ ■ ■ 0/Ja/Vjv + 

5,iVi ■ • ■ ~ ;V! 


(41) 


When computing the full expression given by Eq. (34) for the 
maximal number of E-matrices, one verifies that only one of 
the irreducible components of couples to the Yukawa 

term. Indeed, for N odd, one has k - I and therefore 
only couples to the Yukawa when k = +(-), 

otherwise, for N even, one has k + I and therefore only 
0^i -,«(0,i -/jiv) couples to the Yukawa when k - -(-h). One then 
concludes that the computation of the Yukawa given by Eq. (34) 
can be performed directly using the reducible representation 

Y.b ('I'Ll B c-' E[^.r^2 ■■■ r,„] I'I'/L , ( 42 ) 

without the loss of generality. 

Eor illustrative purpose, the antisymmetric tensor represen¬ 
tations have the following dimension in the case of SO(IO): 

~ 10, 4>f2v ~ 45, ~ 120, (pf2v^cr ~ 210, Cpi^yAa-y ~ 126, 

and^^^2,rr ~ 126. 

One can also express the tensor cpf 2 vA...a- in terms of SU(A^) 
tensors. This can be easily computed by expanding the quan¬ 
tity E^Evr 2 ...Ecr <pfivA...cr in terms of the creation and annihilation 
operators [41] as 


r^iYvYA-Ya-(l>fiyA...a- = b\b\b[...bl(pc,CjCt...c„ 

+ {bib\bl...bl^CiCjcyc„ + perms) 

-H {bibjbl...bl4)crcjCt:.c„ + perms) -H ... 

-H {bibjbk...bn-\bl(l)cycjct...c,yc„ + perms) 
+ {bjb jbjc-‘bii(l)kicjCf;...Ci,) 1 


(43) 


^Representations of dimension with m> N are equivalent to represen¬ 
tations with dimension LL)- 
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where = 4 >... 2 j- + i and (^... 5 ,... = 4 >... 2 j... - i(p... 2 j-\.- 

The new tensors in equation Eq. (43) manifest completely anti¬ 
symmetry, i.e., 

<P - Ci - Cj - — - - 5 (44) 

In Appendix C, we compile all the antisymmetric tensors of 
SO(IO) explicitly written in terms of SU(5) representations. 

2.3. Methods and rules 

We compile in this section the rules that are the basis for 
defining the behaviour of the SO^ library. In order to con¬ 
cretise the use of the rules, let us first take a simple bracket 
containing a set of annihilation and creator operators, 

<^\bibjblbibibl\Q) . (45) 

The tensor differentiates upper and lower indices asso- 
ciated with annihilation and creation operator indices, respec¬ 
tively. This distinction is important to obtain a final expression 
with a index structure consistent with SU(A). The computation 
of the bracket from Eq. (45) relies on the use of the relation 
given in Eqs. (9) and (12). We present two different strategies, 
that we call normal ordering and reverse ordering methods. We 
first discuss the reverse ordering method. 

Reverse ordering method 

In this case, we move according to Eq. (9) either the anni¬ 
hilation operators forward to the right until it cancels with | 0 ), 
bi | 0 )= 0 , or the creation operators backward to the left to cancel 
with (0|, i.e., (0| b\ - 0. Hence, moving the operator bi to the 
right, we have 

< 0 | bibjblbiblbl | 0 > = 

^{Q\bibjbl{6,„-bib)bl\Q) 


= 6l,n <0| bi bj bl bl 10) - 6,n <0| bi bj bl bl 10) . 

We complete the computation of (0| bi bj b^, bl |0) and 
( 0 | bi bj bl bh | 0 ) using the same procedure, obtaining 

< 0 | bi bjblbl\Q)^ 6 jk 6i„ - 6j„ 6 it , 

< 0 | bi bj bl bl^ | 0 ) = 6jk dim - dj,„ da,. 

The final result for (0| bi bj bl bi bl, bl |0), in terms of d’s, is 
then given by 


dlmidjkdm djndik^ dlf,{djkdim djmdif,^. 


(48) 


Normal ordering method 

In the normal ordering method, we use the relations in Eq. (9) 
to rearrange the creation and annihilation operators in such a 
way that all annihilation operators are on the left-handed side 
while all creation operators are on the right-handed side, as 

< 0 | bibjblbiblblm^ 

^{0\bibj{dM-bib^^blbl\0) (49) 

= dki <0| bi bj bi bl 10) - <0| bi bj bl bl bi bl 10) . 


Then we compute (0| bi bj bi bl bi bl |0) and (0| bi bj bl, bl |0) 
by using the relation. 


\bi,bi,...bi,b\bi...b]\0)^ 


1 


(N-ky. 


iii2...ikk+i--iN c.. ... 

^Jk--JlJ\k+\—iN 


(50) 


that holds for SU(A) with k < N. 

Up to now, we did not mention in which SU(A) framework 
we are computing the expression given in Eq. (45). If we choose 
to compute it in SU(5) (i.e., A = 5) the brackets in Eq. (49) take 
the values 


<0| bi bj bi bl |0) = ^ Snmafly , 

< 0 | bi bj b, bl bi bl | 0 ) = i e„mkap , 
and the tensor defined in Eq. (45) becomes then 


(51) 


•y'T - lx' Jjafiy 

^kmn - 


^nmafSy 


enmkafi- 


(52) 


We are ready to summarise all the rules mentioned above. It 
is known that a complete bracket expression is composed by 
creation and annihilation operators, and some other fields. The 
idea is to get rid of all operators by using the relations in Eq. (9) 
and use the terms arising from these calculations to simplify 
the remaining expression (e.g., fields if there are any). In or¬ 
der to have a consistent expression, each element in the bracket 
must obey certain rules, those that we list below for a general 
SO(2A); 

1. The number of creation operators in a bracket expression 
must be equal to the number of annihilation operators. 

2. The number of contiguous creation (or annihilation) oper¬ 
ators must be equal or less than N for SO(2N). 

3. The operators bi inside the bracket expression are written 
on the left-handed side while the operators b] are written 
on the right-handed side, otherwise the result is zero, i.e., 
bi | 0 ) = 0 and < 0 | bl = 0 . 

4. The difference between the number of upper and lower in¬ 
dices in the fields must be zero or multiple of N for any 
SO(2A). 


As already said, instead we move the annihilation operators 
to the right-handed side to cancel at | 0 ), we can move the cre¬ 
ation operators to the left to cancel when reach (0|. However, 
once we choose to move either annihilation or creation opera¬ 
tors, we need to maintain this choice until the end of the com¬ 
putation. 


Within this framework, we shall give the Hermitian conjuga¬ 
tion and the transposition operations on a general vector 


I'T) = 10 ) ^ + b] 10 ) + -bib] 10 ) 




IJ 


( 53 ) 
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[< bra I ket >] 


Braket inl,in2,out: 
setDim(4): 

inl=bra(0,N(A)*e_(i,j)/2, bQ) * b(i)): 
in2=ket{0.Nb(B)*e_(k,l)/2, bt(k)*bt(l)): 

out = ini * in2: 




Figure 1: Scheme of SO^ library exemplifying how to compute 

"-^NaNb <0| bjb.blb] |0> in S0(4). 


User 

Braket/FORM 
BraketOneTerm 


□ List 

elemType / noList 


Figure 2: Hierarchy of the SC^ii library. 

{bra\ ■ free 
free ■ \ket) 

{bra\ + {bra\ 

\ket) + \ket) 
free + free 
free ■ free 


which are 


• {bra\ket) + {bra\ket) 


<'pi = r{o\ + iifi{o\bi + <oi bjbi 

■t ^ ^ {0| bfjiblbjf + , 


(54) 


and 


= ^<01 + <0| bjbi + b,„b,bk 

(55) 

respectively, where we make the usual identification of lower 
indices as tf/ij - (i/r'-') . 


3. SO^, a C++ library 

In this section we present the structure of the SO^ library, 
that can be found in http;//sospin.hepforge.org, com¬ 
ment on the data structure representation and give a list of the 
most important functions to use when writing a program linked 
to the SO^ii library. As mentioned before, the SO^ii code is 
entirely written in C++ and it is based on operations over cre¬ 
ation and annihilation operators. We give in Fig. 1 the pictorial 
scheme to explain how a SO^ii program works. It works like 
this: the building blocks of the code, i.e., the {bra\ and the \ket) 
entities, can be defined either in a main file program (ini and 
in2 in the picture) or included in some include file and called 
in the main file. The implemented operations in the SO^ii li¬ 
brary are the following: 

• {bra\ ■ \ket) 


• {bra\ket) ■ {bm\ket) [This operation is only allowed after 
the evaluation of the expression.] 

Once the expression to evaluate is defined (out in the picture) 
the approach to solve it is the following: the expression to eval¬ 
uate is split in two parts, one with all constants and another one 
containing operators; the operator part will be worked out us- 
i|ig one of the two methods described in Section 2.1 leading to 
an intermediate expression written in terms of e’s or 5’s. Then, 
the constant part and the intermediate expression are joined to¬ 
gether to lead to a semi-final expression. If the starting example 
is simple, the result will be simple. On the other hand, if the 
example is somewhat complex, the expression obtained at this 
stage is rather large and needs extra simplifications for increas¬ 
ing the readability of the final expressions. Hence, in order to 
make the reading of all results as easier as possible, we include 
the possibility to simplify the expression obtained so far with 
the Symbolic Manipulation System FORM [49]. Once the con¬ 
stant part and the intermediate expression are joined together, 
it is created an input file to be run by FORM leading to a more 
simplified expression. The final expression is then read back 
from the FORM output file to the program. Note that, in the 
output FORM file, we first present all partial results and then 
the complete result at the end of the file. 

We have opted to include the FORM program as a tool to do 
the final simplifications if needed. For the sake of curiosity we 
have implemented in FORM all the procedure used in SO^; 
however, the running time measured is far larger than in the 
case of using our library in C++. This shows the performance 
power of our choice for the appropriate data structure, which is 
described in the next section. 
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In Fig. 2 we define how the class structure works. In terms 
of level abstraction the low level implementation of the ba¬ 
sic elements for the expression evaluation are the structures 
elemType and noList. They represent the expression ele¬ 
ments or nodes. The class DList is then build over this abstrac¬ 
tion level and represents a linked connection of nodes. Then 
BraketOneTerm and Braket classes represent complete ex¬ 
pressions that can be evaluated and simplified by FORM and 
by the User at the higher abstraction level. 

3.1. Data structure representation 

In order to manipulate sequences of operators or b], we 
need to find the adequate data structure to store and further 
evaluate such sequences. Such data structures should require 
the following criteria: 

• optimize memory usage - since the sequences can get ex¬ 
tremely long; 

• optimize flexibility of permutations - adjacency in memory 
is not relevant; 

• standardize description of all elements - to ease interpreta¬ 
tion and evaluation, and to reduce memory waste in con¬ 
traction and expansion operations. 

We have adopted the doubly-linked list scheme as the appro¬ 
priate solution to the problem. A doubly-linked list consists 
on the list of connected nodes, which include specific data ob¬ 
jects, such as each node is linked to previous and next nodes 
in the list. This is advantageous since one changes only the 
pointers without modifying the content and their position on 
the memory. The doubly-linked list scheme is implemented by 
a C-n- class named as DList. The full method list is included 
in the Appendix D. 

The different types of elements within a sequence are en¬ 
coded as bit-fields of an integer type (int) with the purpose 
of optimising the memory usage. Thus, each DList-node can 
account for the operators bi and 1?.^ as well as the constants and 
the Kronecker symbol 6, with its indices. In Fig. 3, we illustrate 
the concept of the DList class for a simple sequence. 


the program needs a total of 454.15 MB and 2.80 s to evaluate 
this expression in the delta form and 1.41 MB and 0.00064 s 
to evaluate the same expression to the Levi-Civita tensor form. 
The large amount of memory used to evaluate in the delta form 
is due to the number of terms generated in this way, a total 
of 9! = 362880 terms, while for the evaluation to Levi-Civita 
tensor the result has only one final term. 

3.2. General functions 

In this section we list the most general functions needed to 
write a program using SO^ library; it is divide in three subsec¬ 
tions: generic and building functions as well as specific func¬ 
tions to interface with FORM. 

Generic functions 

• void setDimCint n) 

Sets the group dimension. 

• int getDimO 

Gets the group dimension. 

• void CleanGlobalDecK) 

Cleans all tables with indices and function declarations. 

• void setVerbosity ( Verbosity verb) 

Sets verbosity level; verbosity options: SILENT, SUM¬ 
MARIZE, VERBOSE, DEBUG.VERBOSE. 

• Verbosity get Verbosity () 

Returns current verbosity level. 

Building functions 

• DList b(i)/DList bb(i) 

Declares a operator f>, ; the index in bbO must be enclosed 
in quotation marks or passed as a std:: string type. 

• DList bt(i) /DList bbt(i) 

Declares a operator f>!; the index in bbtO must be en¬ 
closed in quotation marks or passed as a std::string type. 


Computation Performance. In order to give an estimate of 
the performance of the computation, we measured the time 
consumed and the memory used in computing a sequence 
of creation and annihilation operators. We run the test in a 
x64 LINUX machine (Ubuntu) with an Intel(R) Core(TM) i5- 
3317U CPU @ L70GHZ. 

We tested the following expression for SO(IO), 

<01 b,b,b]b,b,A.b]j^Mj’‘MMA 10) - (56) 


the program needs a total of 3.57 MB and 0.0210 s to evaluate 
this expression in the delta form and 1.39 MB and 0.0044 s to 
evaluate the same expression to Levi-Civita tensor form. 
Running the following sequence in SO(18), 


(Oift,,^.1,,.ft,I,;,I,;,6],io>, (5?) 


• DList delta (i,j) 

Declares the dij function. 

• DList identity 

Declares the identity matrix. 

• Braket bra(A, B, C) 

Braket ket(A, B, C) 

Braket braket (A, B, C) 

Braket free (A, B, C) 

The element A corresponds to the global index, B to the 
constant part (e.g. fields) and C to the operators b and b\ 
6 or the identity. The first entry is the sum of the num¬ 
ber of all upper indices (positive counting) and lower in¬ 
dices (negative counting) present in the fields defined in 
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*end 

‘actual 

sign 
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b] 


btj 


bk 

*nxt 

*nxt 

*nxt 



*prv 

*prv 

*prv 




noList noList noList 


Figure 3: Representation of the bracket {0|fo,foj6t|0) in a linked data structure. 


the function. This entry can be set to zero and if so we 
must call first the function unsetSimplifyIndexSumO . 

• void evaluate (bool onlydeltas=true) 

Evaluates expression, if onlydeltas is true then the expres¬ 
sion is evaluated to deltas, if false the expression is evalu¬ 
ated to Levi-Civita tensors with eventual d’s. 

• Braket Bop ( std ;; string startid="i") 

Returns the operator B using generic indices. 

Braket BopIdnumO 

Returns the operator B using numeric indices. 

• void newld (string i) 

Declares a new index. 

• void setSimplifyIndexSum () 
void unsetSimplifyIndexSumO 

Activates/Deactivates internal simplifications based on the 
Braket Index sum. This option is activated by default. 

Specific functions to interface with FORM 

• std;: string Field (A, B, C, D) 

This function is used to declare the field in FORM, where 

A field name; 

B number of upper indices; 

C number of lower indices; 

D field properties: 

- SYM; symmetric field without flavor index; 

- ASYM: antisymmetric field without flavor in¬ 
dex; 

- SYM_WITH_FLAVOR: symmetric field with 
flavor index; 

- ASYM_WITH_FFAVOR: antisymmetric field 
with flavor index. 


Returns field name as it should be written in the con¬ 
stant Braket part. 

The convention to write a field in the constant part of a 
Braket is the following: for each field we assign a name, 
then we write the number of upper indices followed by 
the number of lower indices, then between parentheses we 
add the indices, the first index is always reserved for fla¬ 
vor if applicable, then we write the upper indices by the 
order they appear (left to right) followed by the lower in¬ 
dices (left to right). In the case we have some ambiguity 
concerning the symmetric or antisymmetric nature of the 
indices, we add the s letter for those that are symmetric just 
after the field name, e.g. the field with symmetric i, 
j indices and flavor index b, must be written in the Braket 
constant part as Ms02 (b, i , j) and declared to FORM as 
Field (M, 0, 2, SYM_WITH_FLAVOR). 

• void CallForm (Braket feexp, bool print=true, 
bool all=true, string newidlabel="j ") 

Creates the input file for FORM, run the FORM program 
and returns the result to an output file and/or to the screen. 

• void setFormRenumber () 

Sets "renumber 1;" in FORM input file. This option is used 
to renumber indices in order to allow further simplifica¬ 
tions. However, in large expressions this must be avoided 
since it increases the computational time in FORM. The 
best way to use it is simplify the expression with FORM 
with this option unset, and then send a second time to 
FORM with this option active. By default this option is 
unset. 

void unsetFormRenumberO 

Unsets "renumber 1;" in FORM input file. 

• void setFormIndexSum () 
void unsetFormIndexSumO 
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Sets/Unsets the index sum in input FORM file. The set 
option is activated by default. 

4. Work with SO^ 

In this section, we describe the installation of the SO^ li¬ 
brary and other tools that may be provided for the library to 
make further simplifications. We give in detail instructions how 
to use the SO^ library in an standalone C-i-H- program with a 
very simple example involving the group SO(4) just for illus¬ 
tration. 

4.1. Download and installation 

The SO^ library project is hosted by Hepforge at http: 
//sospin. hepforge. org under a GNU Lesser general public 
license. 

The simplest way to compile the SO^ library is: 

1. ./configure --prefix= 
library_installation_path --with-form= 
F0RM_path 

The user can omit the FORM path declaration and set it 
after using export PATH_T0_F0RM=F0RM_path or put 
the binary FORM file in the folder where the user has its 
project. 

2. make 

3. make install 

4. make doxygen-doc (optional - it generates SO^ library 
documentation) 

Inside the library folder there are several example files, in ad¬ 
dition to the ones shown in this paper, to help with the use of 
the library. The FORM [49] binary files can be downloaded at 
http: //www.nikhef .nl/~form/, after accepting the license 
agreement. They are available for LINUX (32-bits or 64-bits), 
Cygwin (32-bits) and Apple/Intel platforms. All the informa¬ 
tion concerning the installation is written in the README file. 
These procedures were successfully tested in LINUX and Mac 
OS X 10.10 (Yosemite). 

4.2. Writing the first program 

In this section we discuss how to write a first program exam¬ 
ple using the SO^ library properly for SO(4). Although the 
group SO(4) ^ SU(2) x SU(2) is not particularly interesting 
for GUTs, it is invoked to illustrate the use of the library in a 
simpler way. Since SO(4) belongs to the family group SO(2N) 
for N even, the two spinors in which the 4-dimensional space is 
broken are not related by conjugation. The general ket is given 

by, 

Ifi) = |0) M + bl |0) N'^ + y^b]b] |0) M , (58) 

where M,M ~ ^ and A' ~ 2 in SU(2). So, the general ket in 
Eq. (58) can be decomposed as, 

\fi) = -H \fi2) , (59) 


where Fq Ii/'i) = \4f\) and Fq \ 4 i 2 ) - - 1 *^ 2 ), which are written as 

\fix)^\0)M+]^s'jb\b]\0)M, (60) 

\fi2)^bl\0)NK (61) 

Using the rules compiled in Section 2.3 for the transposition 
operations, we obtain the following expressions 

{fi\\^M{0\ + ]^s'"'M{0\b,„bi, (62) 

<^*| =A"<0|Z7„. (63) 

As our first example program, we will address the calcula¬ 

tion of I B \ In order to better understand how the li¬ 
brary works, let us start by showing first how to code | ^ 

ignoring the fields and taking into account only the creation and 
annihilation operators, i.e., 

(<0| + <0|/7„Z7,)(|0> + /7,!/7];|0)). (64) 

To properly write the code for this expression, we first need 

to add the header file for the SO^ library and the sospin 
namespace. After creating the main function we must define 
the group dimension in the beginning, setDim(4), and clean 
up all the memory allocated using CleanGlobalDecl () before 
exiting from the program. The code to solve the problem in 
Eq. (64) is the following 

#include <sospin/son.h> 
using namespace sospin ; 

int main(int argc , char *argv[]){ 
setDim (4) ; 

Braket left = identity; 

left -i-= b(m)*b( 1) ; 

Braket right = identity ; 
right -l-= bt( i )*bt( j ) ; 

Braket res = left * right; 
res. evaluate (); 

res . setON () ; 

s td : : cout«" Re suit : \ n "« res «std : : endl ; 
res . setOFF () ; 

CleanGlobalDecl 0 ; 

exit (0) ; 

1 

The evaluate () function can be used with or without the ar¬ 
guments; true or false. The evaluate ( false ) sets on the results 
written in terms of the Levi-Civita and it is only used in opera¬ 
tions with braket type; evaluate O/evaluate (true) does the 
evaluation to delta functions and it can be used in braO/ket () 
/braket 0 and none types. If we use more than one evaluation 
process, we need to maintain the type of evaluation chosen. 

The functions setONO and setOFFO make possible the 
writing of Local R?= for each term in Braket expressions, in 
addition, each term of the expression is numbered. Note that in 
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the code written above, due to the way how it was declared, the 
terms need to be joined by using the += operation. 

In order to compile and run the program, the user must pass 
to the C++ compiler the path to the SO^ library and include 
folder, as for example (assuming the GCC compiler), 

g++ -03 -I / Sospin_PATH/include -L/Sospin_PATH 
/lib example . cpp -o example -Isospin 

Running the program above, we will get the following result^. 

Local R1 = +1; 

Local R2 = + bt(i) * bt (j); 

3 Local R3 = + b(m) k b( 1) ; 

Local R4 = + d_(m, j) * d_(l,i) 

.3 - d_ (m, i) * d_ (1 , j ) 

- d_(l,i) * bt(j) * b(m) 

+ b(m) * bt(i) K bt(j) * b(l) 

+ d_(l,j) >i< bt(i) * b(m) ; 

This result is not in agreement with the rules given in Sec¬ 
tion 3. Following those rules and looking at the term above we 
see that all the terms containing operators must vanish. The 
reason why these terms appear in the result above is because 
we never declared the type of left, right and res, hence by 
default all these expressions are of type free (operation none). 
In order to properly solve Eq. (64) one needs to setup explicitly 
the type of each expression; we can declare them as 

left . TypeO = bra ; 

2 right . Type 0 = ket ; 

There is no need to declare the variable res because the op¬ 
eration left*right will automatically setup its type based on 
the product, i.e., the resulting type of res is braket. 

A simpler and more complete way to declare the expressions 
in Eq. (64) is the following, 

Braket left = Braket ( identity , bra); 

2 left += Braket(b(m)=i<b( 1) , bra); 

Braket right = Braket ( identity , ket); 

4 right +- Braket ( bt (i)* bt (j), ket); 

where we use the operation += for each contribution in different 
lines or 

Braket left = Braket ( identity , bra) + Braket( 

b(m)=i<b( 1 ) , bra) ; 

2 Braket right = Braket ( identity , ket) + Braket 

( bt (i)K bt (j) , ket ) ; 

where we use the + operation for terms placed in the same line. 

Running the program with the types properly setup, we get 
the expected result: 

Local R1 = + 1; 

2 Local R2 = + d_(m, j) * d_(l,i) 

- d_(m, i) * d _(1 , j ) ; 

If we wish to evaluate the expression in Eq. (64) in such a way 
that the final result appears written in terms of the Levi-Civita 


^Note that, in order to save space, we altered slightly the aspect of the pro¬ 
gram’s output. 


tensors, we need simply set the argument of the evaluate func¬ 
tion to false as res. evaluate ( false ). The result will be given 
by the output 

Local R1 = +1; 

: Local R2 = +e_(m,l)=i<e_(j ,i); 

As one can see, the final result written in terms of the Levi- 
Civita tensor is equal to the one written in terms of deltas but in 
a much more compact form. So, hereinafter we will just present 
the results written in terms of Levi-Civita tensors even though 
both methods are available. 

Let us now discuss how to include in Eq. (64) the operator B 
ofEq. (A. 19), i.e., 

{{0\ + {0\b„,b,)B[\0} + blb]\0}). (65) 

Writing this code is rather simple because the SO^ library 
already have a function to compute the operator B for any group 
S0{2N), BopO, therefore we only need to add this function to 
the code as: 

Braket res = left * Bop() * right; 

There is no need to set the group dimension in Bop 0 function 
because it is already done through the setDimO in the begin¬ 
ning of the code. If one wants to define the operator B by one¬ 
self, without using the predefined function, we can do it just by 
using the free 0 function and the rules given above. The result 
will be 

Local R1 = l/2*e_(il ,i2)*( -e_(i1 ,i2)*e_(j , i ) 

) ; 

2 Local R2 = l/2=He_(il,i2)>i<( +e_(il,tl)=He_(i2,tl 

) ); 

Local R3 = l/2=He_(il ,i2)*( 

4 +d_(il,i2)*e_(m, l)*e_(j,i) 

-d_(il , i )*e_(m, 1)*e_(j ,i2) 

6 +d_(il,j)*e_(m,l)*e_(i,i2) 

); 

8 Local R4 = l/2=He_(il ,i2)>i<( 

+d_(i2,i)>i<e_(m,l)=He_(j ,il) 

10 -d_ (i2 , j ) * e_ (m, 1) * e_ (i , i 1 ) 

); 

12 Local R5 = l/2*e_(il,i2)*( -e_(m, l)*e_(i2,il) 

); 

These results are quite large and clumsy so in order to sim¬ 
plify them we decided to include FORM as a hnal step. The 
inclusion of FORM is purely aesthetics and does not affect the 
computation procedures. In the case we consider it the last re¬ 
sult become 

Local R1 = +e_(i,j); 

2 Local R2 = -e_(l ,m) ; 

After this introduction, we are ready to compute 
(tAj J BI tAii), i.e., 

(m, <0| + (0| b^b,) B (10) Mh + y^b^b] |0) Mb ), (66) 

where a and b are flavor indices. 
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In order to compute this example, we need to add the fields 
M and M. This is done by using the functions bra (), ket (), 
free 0 and braketO defined in Section 3.2. To account for the 
changes due to the inclusion of the fields we need to substitute 
the codes given above by 

Braket left = bra (0 ,M(a), identity ) ; 

2 left +=bra(0 ,Mb(a)*e_(l ,m)/2 ,b(m)=i<b( 1)) ; 

Braket right = ket (0 ,M(b) , identity ); 

4 right +=ket (0 ,Mb(b) =He_(i , j )/2 , bt (i) Kbt (j )) ; 

where the field M is coded as Mb. The output result is 

Local R1 = M(a)*l/2*e_(il ,i2) *Mb(b) *e_(i ,j ) 
/2*(-e_(il ,i2)*e_(j , i)) ; 

2 Local R2 = M(a)*l/2*e_(il, i2)*M(b)*( +e_(il,tl 
)Ke_(i2,tl) ); 

Local R3 = Mb( a)=i<e_(l ,m)/2>i<l/2>i<e_(il ,12) =i<Mb(b ) 
*e_(i>j)/2*( 

4 +d_(il , i2 ) *e_(m, 1)*e_(j ,1) 

-d_(il ,i)*e_(m, l)=i<e_(j ,12) 

6 +d_(il,j)*e_(m, l)*e_(i,i2) 

); 

8 Local R4 = Mb( a)=He_(l,m)/2>i<l/2>i<e_(il ,12) =i<Mb(b) 
*e_(i ,j)/2*( 

+d_(i2,i)*e_(m,l)=i<e_(j ,11) 

10 -d_ (i2 , j ) * e_ (m, 1 ) * e_ (1 , 11 ) 

); 

12 Local R5 = Mb( a)=i<e_(l ,m)/2>i<l/2>i<e_(il ,12) *M(b) 

*( -e_ (m, 1) >i< e_ (i2 , i 1 ) ); 

To make the final simplification we use the FORM program; 
before we call it to simplify our expression, we first need to 
declare explicitly all fields as well as all indices appearing only 
in the constant part, i.e., the indices a and b in this example. 
Therefore, we need to add the following code. 

Field (M, 0, 0, ASYM_WrrH_FLAVOR) ; 

2 Field (Mb, 0, 0, ASYM_WrrH_FLAVOR) ; 
newld("a"); newld("b"); 

Once the fields M and M carry flavor, we need to set them as 
ASYM_WITH_FLAV0R, for more details see Section 3.2. To call 
the FORM program to simplify our expression we only need to 
write 

CallForm ( res , false , true, "j"); 

Note that the function callFormO sets setOFFO for the ex¬ 
pression Braket. If the user have declared setONO previ¬ 
ously, the user must set setON () again for that expression after 
callFormO . For a more detailed description about this func¬ 
tion please see Section 3.2. 

Running the program above we obtain the following result. 

Local R1 = +M(jl) * Mb(j2); 

2 Local R2 = - Mb(jl) * M(j2); 

Note that jl and j2 are not arguments but indices; they corre¬ 
spond to the flavor indices a and b. This is so because by default 
we consider that the indices are summed and hence they are re¬ 
named. If we want to avoid summed indices we must add the 
function unsetFormIndexSum () and then the result would be 


Local R1 = -H M(a) * Mb(b) ; 

Local R2 = - Mb(a) * M(b) ; 

For the sake of completeness, we give below the complete 
code to compute | B | in SO(4). 

#include <sospin/son.h> 
using namespace sospin ; 

int main(int argc , char *argv[]){ 

setDim (4) ; 

Braket left = bra (0 ,M(a), identity ); 
left -H=bra(0 ,Mb(a)*e_(l ,m)/2 ,b(m)*b( 1)) ; 

Braket right = ket (0 ,M(b) , identity ); 
right +=ket(0 ,Mb(b ) >i<e_(i,j)/2,bt(i)>i<bt(j)); 

Braket res = left * Bop() * right; 
res. evaluate(false) ; 

Field (M, 0, 0, ASYM_WrrH_FLAVOR) ; 

Field (Mb, 0, 0, ASYM_Wrm_FLAVOR) ; 
newld("a"); newld("b"); 

unsetFormIndexSum 0 ; 

CallForm ( res , false , true, "j"); 

res . setON () ; 

s td : : cout «" Re suit : \ n "« res «std : ; endl ; 
res . setOFF () ; 

CleanGlobalDecl 0 ; 

exit(0); 

1 

where the result is obviously Ma Mb - Ma Mb ■ In order to com¬ 
pute I B I tl/ 2 h), ((^2 I BI lAi) and | B \ 1 /^ 2 *) we just need 
to define and \ 4 f 2 b), and substitute it in the code above. 
Using Eq. (63) we define {4'2i\ 

Braket left = bra(l ,N(a,n) ,b(n) ) ; 
and using Eq. (61) we define \ 4 i 2 b) as 
Braket right = ket (1 ,N(b ,k), bt (k)); 

The results are: 

5. Examples 

In this section we present two more complex examples to 
better explain how to use SO^ library. We give one example 
in SO(4) with higher dimensional terms and one example in the 
context of SO(IO) models. 
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16 

18 
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26 

28 
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Table 1: The 90^ code to compute \ BY^ | \l/ 2 b^ {^\c I I in SO(4) and the corresponding result. In this example we have used the operator B, BopO, 
that uses generic indices, i and k. 


#include <sospin/son.h> 
using namespace sospin ; 

4 int main(int argc , char >i<argv[]){ 

6 setDim (4) ; 

8 Braket LI, R1, L2, R2; 

Braket inlE , inlO , in2E , in20 , res; 

10 

LI = bra(0 ,M(a), identity ); 
i: Ll+= bra(0 ,Mb(a)*e_(i , j )/2 ,b( j )*b( i )) ; 
R1 = ket(0,N10(b,k) ,bt(k)); 

14 

L2 = bra(0 ,M(c), identity ); 

16 L2+= bra(0,Mb(c)=He_(l ,m)/2,b(m)=i<b(l)); 

R2 = ket(0,N10(d,o) ,bt(o)); 

18 

Field (M, 0, 0, ASYM_Wrra_ELAVOR) ; 

20 Field (Mb, 0, 0, ASYM_Wrm_FLAVOR) ; 
Field (N, 1, 0, ASYM_WrrH_FLAVOR) ; 


newld ( "a" ) ; 

newld ( 

"b' 

') 




24 newld ( " c " ) ; 

newld ( 

"d' 

') 




26 inlE = LI * 

Bop( " i 

") 

* 

G( true , 

"j") 

=8 Rl: 

in2E = L2 >8 

Bop( "k 

") 

* 

G( true , 

"j") 

=8 R2: 

28 inlO = LI >8 

Bop( " i 

") 

* 

G( false , 

"j") 

=8 Rl: 

in20 = L2 * 

Bop( "k 

") 

* 

G( false , 

"j") 

=8 R2: 


inlE. evaluate () ; 

32 in2E . evaluate () ; 

inlO . evaluate () ; 

34 in20 . evaluate () ; 

36 res = inlE * in2E + inlO * in20 ; 

38 unsetFormIndexSum 0 ; 

CallForm ( res , false , true, "i"); 

40 

res . setON () ; 

42 std : : cout « "Output result :\n" « res « std::endl; 

44 CleanGlobalDecl 0 ; 

46 exit(O); 

) 


Output result: 

Local R1 = +2*M( a) =i<N10(b , i ) =i<Mb( c) =hN 10(d , j ) =8 e_ (i , j ) ; 
Local R2 = -2KMb(a)*N10(b, i)*M(c)=8N10(d, j )=i<e_(i , j ) ; 
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5.1. S0(4) 

Let us compute the following higher dimensional interaction 
term in SO(4); 

^ Z I I (^ 1 -^ I I ^ 2 ^) ’ 

H=\ 

where Yah, Ycd are Yukawa matrices, a, b, c, d are flavor indices 
and A is some high energy scale. The above equation can be 
expanded by rewriting the T^ in terms of the operators bi and b^. 
as given in Eq. (10) 

2 

I ^ I {flc I B{bi + b\) I ilf2d) 

i=i ^ (69) 

- I B{bi - b]) I ilf2t) {ric I B(bi - bj) I l/f2d)) ■ 

The code to compute Eq. (68) is given in Table 1 and the 
corresponding T^ is given by taking into account the parity as 

Braket G(bool even, string startid){ 

2 if (even ) ( 

Braket G_even = bb(startid +"1"); 

4 G_even += bbt( startid +"1"); 
return G_even; 

^ } 

Braket G_odd = bb( startid +"2"); 
s G_odd -= bbt( startid +"2"); 

string constpart = 

10 G_odd = G_odd * constpart; 
return G_odd; 

■2 1 

Note that it is important to carefully define different indices 
among bracket expressions in order to avoid repeated indices, 
which could lead to a meaningless result. The result of the 
above code is then given as 

2^^ (m„A' Sij. (70) 

5.2. SO(IO) 

In this subsection, we give the example 16fl16i120// of 
SO(IO) computed using the SO^ library. In a first step we 
present the program by defining all quantities while in a second 
step we rewrite it using only the specific SO(IO) functions al¬ 
ready included in the library. The reducible 32 representation 
of SO(IO) is given by 

I'D = |0> ^ + b] |0> If/ + ^-b]b] |0> + ^s‘^'“"'blb]bl |0> 

+ |0) + blblblblbl |0) ^ . 

(71) 

As already mentioned, in SO(IO) the fermionic particles are 
usually assigned to the 16 dimensional representation which 
corresponds to the semi-spinor '!'+ while 16 s The spinor 


SO(IO) 

SU(5) 


10) 

1 

M 

i.’lo> 

5 

M‘ 

b]bl 10) 

10 

M'j 

b]blbl\Q) 

To 

Mij 

b]b[b]bl 10) 

5 

Mi 

b\blblblbl 10) 

1 

M 


dim = 2^ = 32 


Table 2: The SO(IO) states in terms of SU(5) fields. 


representations are schematically given in Table 2 where |'1'+) 
and !'!'_) are given by 

I'PD = 10) M + ^-b]bl 10) |o) m . ^ (72) 

and 

I'PD = bl 10) M' + ^s‘^'^"^blb]bl 10) Mij +blblblblbl 10) M, 

(73) 

while the transpose of I'T), represented as {'I'*| are given by, 

{'i'll ^M{0\ + <0| bjbi + (0| bibkbjbi , (74) 

and 

('Dl = M‘ <0| bi + Mij <0| b„Abk + M (0| ^ 75 ^ 74 ^ 3 ^^!, 

(75) 

where the flavor index was omitted. 

The Yukawa term Yah 16a 16;, 1 2Qh is written as 

7? E^DEp !'!'+,) (b^vp. (76) 

In order to compute this expression using the SO^ library, 
we write our main program as done before for SO(4) examples. 
The code to compute this example is given in Table 3. We start 
by including the son.h file, sospin/son.h in line 1, then we 
set the group dimension with the function setDim(lO) (line 5). 
The expression for |'I'+) is given in Eq. (72) and is declared in 
line 11; ('P*| is given in Eq. (74) and declared in line 7 while 
the flavor indices a and b are declared in line 15 in Table 3. 

Making use of the basic theorem [42] pointed out in Eq. (43), 
one can write the action of E^EyEp over the 120-dimensional 
Higgs held (bpyp which is given in Eq. (C.8) of Appendix C and 
defined in line 20 in Table 3. The charge conjugation operator, 
B, is in this code taken as the internal function BopO (defined 
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38 


Table 3: The SO^ code to compute 16a 16i, 120//, i.e., ^ ^ r^ryl^p I'P+i) ^/ivp in SO(IO) and the correspondig result. 


#include <sospin/son.h> 
using namespace sospin ; 

int main(int argc , char >i<argv[]){ 
setDim (10); 

Braket psipbra = bra(0, M(a), identity); 
psipbra += bra(2, 1/2*M20( a , o , p) ,-b (o) *b (p)) ; 

psipbra += bra(4,l/24*e_(o,p,q,r,s) =HMb01 (a,o) ,b(p)=Hb(q)=i<b(r)=i<b(s)) ; 

Braket psipket = ket (0 ,M(b), identity); 

psipket += ket(2,l/2 * M20(b,j ,k), bt(j) * bt(k)); 

psipket += ket(4, 1/24*e_(j ,k, 1 ,m,n)*Mb01(b,j) , bt(k)*bt(1)*bt(m)*bt(n)); 

newld("a"); newld("b"); 

Field (M, 0, 0, ASYM_WrrH_FLAVOR) ; 

Field (M, 2, 0, ASYM_WrrH_FLAVOR) ; 

Field (Mb, 0, 1, ASYM_Wrm_FLAVOR) ; 

Braket gamma = free (-3 ,1/6* ( e_ (rl , r2 , r3 , r4 , r5 ) *H20( r4 , r5 ) / sqrt (3) ) ,b ( r 1 ) *b ( r2 ) *b ( r3 )) ; 
gamma += free(3 ,l/6*(e_(rl ,r2 ,r3 ,r4 ,r5) *H02( r4 ,r5)/sqrt(3)),bt(rl)*bt(r2)*bt(r3)); 
gamma += free(-l ,(2*H12(rl , r2 , r3 )+d_(rl , r2)*H01(r3 )-d_(rl , r3 )*H01(r2 ))/(2* sqrt (3) ) , bt (rl )*b(r2 
)=Hb(r3)) ; 

gamma += free (1 ,(2*H21(rl , r2 , r3 )+d_(rl , r3 ) *H10( r2 )-d_ (r2 , r3 )*H10(rl ) )/(2 * s qr t (3) ) , bt (rl )*bt (r2 
)*b(r3)); 

gamma += free ( -1 ,2*H01 ( r 1 ) / sqrt (3) ,-b ( r 1 )) ; 
gamma += fr ee (1 ,2* H10( r 1 ) / s qrt (3) , bt(rl)); 
newld("r4"); newld("r5"); 


Field (H, 

1 , 

0, 

ASYM) ; 

Field (H, 

0, 

1 , 

ASYM) 

Field (H, 

0, 

2, 

ASYM) ; 

Field (H, 

2, 

0, 

ASYM) 

Field (H, 

2, 

1 , 

ASYM) ; 

Field (H, 

1 , 

2, 

ASYM) 


Braket exp = psipbra * Bop() * gamma * psipket; 
exp. evaluate (); 

CallForm (exp , true , true); 

CleanGlobalDecl 0 ; 

exit(0); 

) 


Output result: 

R = + 2*M(jl )*H10(j2)*Mb01(j3 , j2 ) * sqrt (1/3 ) * i_ 

- M(jl)*H02(j2 ,j3)*M20(j4 ,j2,j3)*sqrt(l/3)*i_ 

+ M20(jl ,j2 , j3 )*H01(j3 )*Mb01(j4 ,j2)*sqrt(l/3)*i_ 

+ M20(jl ,j2 ,j3)*H02(j2 , j 3 ) *M( j 4 ) * s qr t (1/3) * i_ 

+ l/4*M20(jl ,j2 ,j3)*H21(j4 ,j5 ,j6)*M20(j7 , j6 , j 8 ) * sqrt (1/3 ) *e_(j2 , j3 , j4 , j5 , j 8 ) * i_ 

- l/4*M20(jl ,j2 ,j3)*H21(j4 ,j5 ,j6)*M20(j7 , j8 , j6 ) * sqrt (1/3) *e_(j2 , j3 , j4 , j5 , j 8 ) * i_ 

- M20(jl ,j2 ,j3)*H12(j4 ,j2 ,j3)*Mb01(j5 , j4 )* sqrt (1/3 )* i_ 

- 2*Mb01(jl ,j2)*H10(j2)*M(j3)*sqrt(l/3)*i_ 

- l/2*Mb01(jl , j2)*H01(j3 )*M20(j4 ,j2,j3)*sqrt(l/3)*i_ 

+ l/2*Mb01(jl , j2)*H01(j3 )*M20(j4 ,j3,j2)*sqrt(l/3)*i_ 

+ 2*Mb01(jl ,j2)*H20(j2 ,j3)*Mb01(j4 , j3 ) * sqrt (1/3 ) * i_ 

+ Mb01(jl ,j2)*H12(j2 ,j3 ,j4)*M20(j5 , j 3 , j 4 ) * s qrt (1/3 ) * i_ ; 
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in Section 3.2). In line 26 we declare the indices that appear 
in gamma but are not defined automatically (i.e., the ones that 
appear only in fields, r4 and r5 in this case). The definition of 
fields in FORM language is done in lines 16 and 28. In line 32 
the expression is evaluated, in line 34 we call the FORM to 
perform some final simplifications and in line 36 we clean tem¬ 
porary allocated memory with CleanGlobalDecl ( ). 

Note that Yab carries flavor dependency, but it is not written 
explicitly in the code, and hence the final result presented in 
Table 3 can be simplified by performing the symmetrisation of 
the fields using 

A,- Bj = ^(A,- Bj + Aj Bd + ^(A,- Bj - Aj B,)- (77) 

If we work out the final result of Table 3 by using the expression 
above, we obtain 

i ^ y,, (2 MaMb H‘ + m‘J Mb Hij 

+ Mia Mjb H‘j - m‘J Mb Hj (78) 

+ Mia Mf H), - i Sijbbn m‘J MZ' , 

, ^ab ~ ^ba 

where Y ,= -. 

ab 2 

Now we write down the specific function of SO(IO) that are 
implemented in SO^ library; these functions are defined in 
tools/solO. h file. 

• Braket GammaH(int n) 

Gives the action of the n F-matrices acting over the Higgs 
field, as summarised in Appendix C; n is the number of 
F-matrices (runs from 0 to 5). 

• Braket psi_16p(0PMode mode, string id) 

Braket psi_16p(0PMode mode) 

These functions give (i/'* |, defined in Eq. (74), if the mode 
is bra and \ip+), defined in Eq. (72), if we select the ket 
mode. 

• Braket psi_16m(0PMode mode, string id) 

Braket psi_16m(0PMode mode) 

These functions give (i/'ll, defined in Eq. (75), if the mode 
is bra and defined in Eq. (73) if the selected mode is 
ket. 

• Braket psi_144p(0PMode mode) 

Braket psi_144m(0PMode mode) 

Enactions for 144+ and 144_, defined for modes bra and 
ket. The complete expressions are given in Ref [43]. 

In what follows, we rewrite the code presented in Table 3 just 
making use of the included functions described above. 

#include <sospin/son.h> 

#include <sospin/tools/solO.h> 


using namespace sospin ; 
int main(int argc , char *argv[]){ 

setDim (10); 

Braket res = psi_16p(bra) * Bop() * 

GammaH(3) * psi_ 16p ( ket ); 

res. evaluate (); 

CallForm ( res , true, true); 

CleanGlobalDecl (); 

exit(0); 

1 

6. Conclusions 

In this paper we have presented the SO^ library, pro¬ 
vided under the terms of the GNU Lesser General Public Li¬ 
cense as published by the Eree Software Eoundation. This 
is a C-i-H- tool whose main goal is to decompose Yukawa in¬ 
teractions, invariant under SO{2N), in terms of SU(A) fields. 
The library project is hosted by the Hepforge website http: 
//sospin.hepforge.org. 

This library relies on the oscillator expansion formalism that 
consists in expressing the SO(2N) spinor representations in 
terms of creation operators, b] , of a Grassmann algebra, act¬ 
ing on a vacuum state-vector. The SO^ code simulates the 
non-commutativity of the operators and their products via the 
implementation of doubly-linked-list data structures. Such type 
of structures are the ideal method to deal with the usage of long 
chains of products of operators b] and fi,. In this type of im¬ 
plementation, the sequences are linked through nodes that con¬ 
tains information about the previous and the next nodes of the 
sequence; this connections are called links. Moreover, the data 
storage in the memory does not need to be adjacent, this is one 
of the reason why the doubly-linked-lists led to high perfor¬ 
mances in our tests. 

In order to understand the manipulations that SO^ need to 
perform, we reviewed in detail the oscillator expansion formal¬ 
ism. Then, we applied the method for decomposing the SO(2N) 
Yukawa terms with respect to SU(A) interactions. The general 
structure of SO^ library was presented by listing the generic 
and devoted SO(IO) functions. After explaining the installa¬ 
tion and the writing of the first program, we showed in detail 
the usage of SO^ through complete examples in both SO(4) 
and SO(IO) frameworks. Additionally, we provided an higher 
dimensional field-operator example computed in the context of 
SO(4) to illustrate how such a term can be processed with this 
library. Einally, we described the functions available in SO^ 
that were made to simplify the writing of spinors and their in¬ 
teractions specifically for SO(IO) models. The code includes 
also functions to deal directly with the 144 and 144 represen¬ 
tations of SO(IO) and one can compute other quantities beyond 
the results already computed in Ref. [43]. 

We are planning to enhance the use of the memory in our li¬ 
brary by implementing new forms of simplifications. We intend 
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to make the simplifications of the final expressions independent 
of external programs in order to reach more performance. Our 
future plans also include the extension of the library with spe¬ 
cific functions, which are now only implemented for SO(IO), 
automatised for a generic SO(2A^). Furthermore, although the 
SO^ library was projected to cover the groups SO{2N), it is 
easily adapted to groups SO{2N H-1) or other algebraic systems 
that are described by creation and annihilation operators. 
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Appendix A. Clifford algebrae and SO(2A) 

In what follows we cover the Clifford algebra discussion 
given in Ref. [50] and we adapt it for the specihc algebra ob¬ 
tained from SO{2N). In general, we dehne the Clifford algebra 
as 


The exponent is obtained through the general formula for the 
sum of p elements of an arithmetic progression mi, ..., Uj„ i.e., 
S p - p '2 ^^^. The hrst matrix of the right-handed side 
of Eq. (A.6) has to anticommute p - I times, the process ends 
with El, which does not need to move. One may then dehne 
the following projectors in analogy to the chiral projectors for 
fermion helds: 



Pr = 


1 + To 
2 


(A.7) 


They are indeed projectors, since 

Pi + Pr = 1, Pi = Pl, pI^ Pr, 
PlPr = PrPl = 0. 


In order to deduce the main results of this appendix, construct 
the set Cd of all products of T-matrices of the Clifford algebra. 


Cd - {±1, iTm., ±PmXm2, •■■>±ro), 


(A.9) 


which is a hnite group. Thus, we can enumerate all irreducible 
representations of the group Cd compatible with the algebra. In 
particular, the unidimensional representations of Co cannot sat¬ 
isfy the anticommutation of Eq. (A.l). The order of the group 
Co is given by 


IN 


Ord(Co)= 22P 
/=0 ' * 


= 2 


2V+1 


(A. 10) 


{r^^Ev} = 2r]p6py, (A.l) 

where rjp is a real constant with \rip\ - 1. From Eq. (A.l) one 
getsE^ = rjpl and 

Tl^r^pYp. (A.2) 

For the spinor representation of SO{2N), the corresponding 
Clifford algebra has rjp - I for any index p, i.e., 

[TpSv] ^2Spy. (A.3) 

Moreover, one may construct a new independent matrix Eq, de- 
hned as 

Eo = i'^EiEz-'-Ez^, (A.4) 


The hnite cardinality of the group Co ensures that there exists 
a basis where all matrices E^ and their products can be taken 
unitary. From Eq. (A.2), one has in this basis E^ = E^ and 

pT _ r* 

For any hnite group, the number of irreducible representa¬ 
tions is equal to the number of classes. Thus, the enumeration 
of classes, ric, of Co is: 


[1], (-1), l+Ei), l+Ez), ■■■,{±E2 ^), 
, {+E 1 E 2 ), ■■■ l+Eo), 


which implies that 



(A. 11) 


(A. 12) 


which anticommutes with all E-matrices. This dehnition im¬ 
plies that Eq = Eo, since 

rj = (-*rr2^^ ■ • ■ rr^ = /'^(-i)'^E2« ■ • ■ Ei 

•N / 1 \N(2.N—1)+Nt^ Y' 


where one has used the following relation for non-repeating p 
E-matrices with p < 2N: 

rm^-'-r^jEm, = (-i)^e,„,E;„ 2 •■•rmp ■ (a.6) 


We conclude that the hnite group Co has 1 - 1 - 2^^ irreducible 
and nonequivalent representations. In order to determine how 
many unidimensional representations of Co exist, we recall 
that in any hnite group, the number of unidimensional rep¬ 
resentations is given by the ratio of number of the elements 
of the group and the elements of its commutator subgroup. 
The commutator subgroup of Co has two elements, namely 
[Co, Col = (1,-1). Thus, the number of unidimensional rep¬ 
resentations is #Co/#[Co, Co] = 2^^. We then obtain that there 
is only one non-unidimensional representation for the group 
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Co- The non-trivial representation can be obtained through the The product B*B then commutes with all elements of the spinor 
relation representation U{to). The first Schur lemma obliges that 


1^ + 1^ + ■ ■ ■ + 1^, +n2 = 2^^^^, (A.13) 

2™ 

which implies that « = 2^. It is clear that all representations 
of the Clifford algebra are also representations of Co and there¬ 
fore the Clifford algebra has only one irreducible representation 
with dimension 2^. It is interesting to verify that from a given 
irreducible representation of the Clifford algebra T^, the matri¬ 
ces r*, r^, and -T^ form also an irreducible representation of 
the algebra since 

{r;,r:} = {r;^,rj} = {-r;i,-rj} = 2^- (a.i4) 

Due to the fact there is only one irreducible representation of the 
Clifford algebra, these three representations must be necessarily 
equivalents. Thus, there must exist matrices B and C such that 

= (A.isa) 

= (A.isb) 

c^‘r;ic = -r„. (a.isc) 


B*B = e'l, (A.22) 

with e' - +1, since |B| = 1 and taking the conjugation 
ofEq. (A.22) 

e'* = BB* = BB*BB ^ = e'. (A.23) 

The relation between e' and e can be established by observing 
that (T^)* = (rp'’', which implies 

[b'B, T^] = 0, (A.24) 

The combination B^B commutes with all elements of the spinor 
representation U{a>). The hrst Schur lemma leads to B^B - A1 
with A - 1. This can be obtained since |B| = 1 (implying 
|/l| = 1) and for any vector |v) with norm 1 one has 

4 = = |B|v>|^ > 0. (A.25) 

We then obtain that B is unitary and as a consequence that C is 
also unitary. Moreover, using the unitarity of B, 


Note that the operator B is the same for both Eqs. (A. 15a) 
and (A. 15b) since E^ = E* and the operator C can be written as 


C - BEo, 


(A.16) 


1 = BT5* ^ ee'BB ', (A.26) 

one concludes that e' - e. 

One may raise the question how to relate the sign of e with 
the dimension of SO(2A). We start by noting that 


apart from an overall complex factor. This derivation ensures 
the existence of the operator B, which was used in Eq. (26). 

We end this appendix by showing that B and C can be 
taken unitary and they are either antisymmetric or symmetric 
matrices, depending of the value of N. Taking the fact that 
(EpT = and using Eq. (A. 15a) one concludes 

[B-^B^yW^{B-^B^)^r^, (A.17) 

which implies that B^^B'^ commutes with all elements of the 

spinor representation U{oS) and therefore by the first Schur 
lemma one has 

B-^B^^e\. (A.18) 

Without the loss of generality, we can refactor the matrix B to 
have |B| = 1. This choice implies that |e| = 1 and 

B^^eB. (A. 19) 

The constant e is in fact real, since transposing Eq. (A.18) one 
has 


rj = BTqB-^ = CEqC^', (A.27) 

and one deduces, using Eq. (A.16), that 


CT = eC. 


(A.28) 


If one takes into account Eq. (A.6), one is able to write the 
following relation 

/ \T b(p+1) 

= e(-i)—( a.29) 

for p < 2N. The above relation implies that the matrices 
CE^jE^j ••■f’/ip are either symmetric or antisymmetric matri¬ 
ces. On the other hand, the set {CE^jE^j • ■ ■ E^^) forms a basis 
of the vector space of all 2^ x 2^ complex matrices, since all 
order product of distinct E-matrices are 2^ linearly independent 
matrices. Thus, the number of independent antisymmetric ma¬ 
trices is given by 


IN , 

Z 


p=0 


2N 

P 


(A.30) 


e = B{B-^y = B(BT)-i = e*. (A.20) 

Thus, one has e = ±1, where the sign will depend only on N. 
Instead, if we use the relation (Ep* = E^, we get 

[b*B, E^]=0. (A.21) 


Eor any 2^ x 2^ complex matrices, the number of independent 
antisymmetric matrices is just given by 2^(2^ - l)/2, therefore 
one has 


2N 


p =0 


Z2l>-'<-»‘ 


2N' 

P 


2^(2^ - 1) 


(A.31) 


18 






The above equation gives a closed relation between e and N. 
After some algebraic simplibcations, one gets 


2N 


e - 


1 , (p-2)(p-i) /2A 

ov+i ^ 

/;=() 


Taking into account the following relation 

<p-2)(p-l) 1 


(-!)■ 


(l + i)F + (1 -OHT 


(A.32) 


(A.33) 


that can be verified by mathematical induction, one derives 


The fact that two distinct pair of T-matrices anticommutes, it 
guarantees that the operators bj and b\, thus constructed, satisfy 
fully the Grasmann algebra, i.e., 

[bj,bl]^6j,\, [bj,b,] = [b],bl} = 0. (B. 8 ) 

We have shown that for any Clifford algebra one can write a set 
of creation and annihilation operators and therefore the SO^ 
library can be easily adapt for those cases. 

Appendix C. SO(IO) compendium 


V2cos ^( 2 A+ 1 ). (A.34) 

This equation has the period N ^ N + 4, and implies that e = 1 
for A = 3, 4 (mod 4), while e = -1 for A = 1, 2 (mod 4). For 
instance, in the case of SO(IO), one has B'^ - -B and = -C. 

Appendix B. Clifford vs. Grassmann algebrae 


As a matter of completeness, in this section we compile 
some SO(IO) details and functions not written in the previ¬ 
ous chapters. In particular we summarise [41, 42] the action 
of F-matrices on the Higgs fields, (^, as well as the action of the 
operator B on | and (i/^H . 

The generic expression for the operator B is given in Eq. (28) 
and in what follows we list the action of the operator B on (i/'* | 
and \ given in Eqs. (72) to (75), respectively. 


In this appendix, we show the existence of a one-to-one cor¬ 
respondence between Clifford and Grassmann algebrae, apart 
from an overall complex phases. In order to demonstrate this 
result, let us take an arbitrary F-matrix pair F^ and Fy of an 
abritrary Clifford algebra given in Eq. (A.l). Without the loss 
of generality, we shall take two consecutive gamma matrices, 

{r 2 ,-i,F 24 = o, (b.i) 

We identify then two possibilities 772,-1 = 7727 and 7727-1 = -'Hij- 
It is then straightforward to verify that there exist a linear com¬ 
bination bj - Q'F 27 -i + l3T2j such that they generate a Grass¬ 
mann algebra, 

[bj,b\]^\, b]^0. (B.2) 

Thus, when ij = 7727-1 = 7727 one has 

= ^ 0 T 27 -i +F 27 ), (B.3) 


<T':| B = - 7 <0| b„ - <0| bkbib,„ 

- 7 l/f (0| hih2^3^4^5 , 

and 

(T'l| B = (0| bjbkb,b„ + '^^^j <0| b„bj + 7><0| . 

(C.2) 

One sees from the equations above that the action of B on the 
bras { 1 / 7 * I and ( 1 / 7 * | gives expressions which are consistent with 
the SU(5) convention for upper and lower indices. 

We compile below the action of F-matrices on Higgs repre¬ 
sentations of different dimensions, which follows from Eq. (43) 
and Refs. [41, 42]. Eor the 10-dim Higgs representation, (p^, 
one has, 

= bl (f>c„ + b„ (pc,,, (C.3) 


which then implies 



with the following normalisation coming from the SU(5) nota¬ 
tion 

(B.4) 

(pc, = V2A", (Pc, = V2//„. (C.4) 


where we have used the result that F], = rja^a . The two 
Eqs. (B.3) and (B.4) can be easily inverted so that the r 27 -i 
and F 27 are written in terms of bj and hj as 


Eor the 45-dim representation [42], (p^y, one has, 

= bibj(pc-c. + b]b\pc,cj + 2blbj(pc,cj - 4>c„c„, 


r 27 -i = -i(bj - Tjb]), F 27 = bj + rjb]. 


(B.5) where 


In the case of having 7727-1 = - 7727 , one has instead 

bj = i(F27-i +F27), b ] = ^ (F27-1 - F27) , 
or the inverted system of equations 
r27-i = bj - rjb ], F27 = bj + rjb ]. 


4’c„c, — h , pacj — h‘j + \S'jh , 
4^ciCj ~ b ^, ^CjCj ~ bjj , 

(B.6) 

with the normalisation 

/7=VI0A, 

(B.7) hij^yllHij, /7' = V2/B. 


(C.5) 


(C. 6 ) 


(C.7) 
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For the 120-dim Higgs field [42], one has, 




+ 3(bjbjbk(l)ciCjCk + blb'^jbk(l)ciCjCk) 

+ i3b i(pCnC„Ci 3b 1 (pCiiCnCi) ■> 

(C. 8 ) 

where 


<f>CiCjct = - d^/z ) , 

(C.9) 

4>CiCjCi, - bji^- - (bjhk - , 

(C.IO) 


(C.ll) 

^CjCjCif ~ ^ijklmb , 

(C.12) 

0c„c,jCi — b , 

(C.13) 

0c„c„Ci — ~bi , 

(C.14) 

with the following normalisation 


/z' - — h- - — H 

n V 3 ’ ' V 3 ' ’ 

h'j - _L zyU k - — H 

/z'-' = — h' = — H‘ 

"k V3 k ’ '^Jk V3 jk ■ 

(C.15) 


Note that we have corrected the signs of the tensors 
inEqs. (CIO) and (C.14). 

For the 210-dim representation [42], (p^vpAi one has, 

F^FyFpF,i0^yp,, ^Ab\b\blbi(l)c,cjctc, + ‘^b\bjbkbi^cfcfctc, 

+ b]b^jblb\(f>c,cjctc, + bibjbkbicpcfcjctc, 

- 6b]b^jCl)ciCjc„,c„, +6bibj(l)ciCjc,„c„ (C.16) 

+ 3 (f>c„C„,C„C„ ~ ^'3‘b^lbj4>c,CjC,„C,„ 

+ 6b]b^jbkbi(f>c,cjCtc,, 

where 



- — h 

C„jC„C,j '*'5 


(pCiCjC,„C,„ 

= k)*\Sih. 


(pc 

'-iCjCkCi — 2^^Uklmb , 

(pCiCjCkCi 

- ^ ijklm y 

“ 24^ 


(Pci! 

- - h'j 

CjCifjCifj — ? 


(pCiCjC„,C„, 

— bij. 


(pCiCjCkCi — b^i + 

+ 

-is; 

1 


(C.17) 


1 

+ - 


didi) h, 
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k 1J ’ 



(pc, 

- -/z'^-t- 

iCjCkCi ^ 




4^CiC 

- - - w + 
jCkC} - rtjkl ^ 

2 {^‘ibjk 

- 6[hji + 

^jbjk) J 


with the normalisation 




h - 

-- ^-4h, 

V3 ’ 

h‘ = 8 

i V6//' , 

hi = 8 ^/6Hi , 


h'j 

= V2//9, 

bij — 

V2H,i, 


(C.18) 

h‘^‘ 

"i 

k _ V2 TTijk 

b‘u = 
jkl 

Vijji 

•cB jkl ^ 

II 



Finally, for the irreducible representations 126 (cfipvApcr) and 
126(0p yApcr), as it was observed in Section 2.2, one needs only 
to take into account the reducible 252 representation, since only 
one of its irreducible components survives [41], 


r/iryF^FpFp- Apy.|p(y — /Z 

_i_ Ji + 13 b^-W 

- 20 b]blb„h‘ + 5b]b%blb„,h‘ 

+ 15 bihi - 20 blbnbihi + 5 blb„b]j}j}ihi) 

+ 10 /!'■''= - b]b]blb%h‘^'‘ 

+ bibjbkhijk - blb„bibjbkhijk) 

+ (60 b]b]bkh‘j - 30 b]b]bkblbnh‘i 
+ 60 b'lbjbkh), - 30blb„b]bjbkh),) 

+ (5 blb]blb\bJl^‘ + 5 b]bjbkb,bJi‘j,J 
+ (I0b]b]blbibj/^ + 10 b]b]bkbib^h‘ij, 
where 


h = 


VT5 


H, 


4 V2 . 

h' = , 


V5 




V2 

Vis' 


J jklm 




jklm ^ 


*■ jklmn ^ ( 5 ") 5 


2 V 2 


fijk _ 
Im 


VI5 




ijk 

Im 


(C.19) 


(C.20) 


The field denotes the 15 representation of SU(5), which is 
a symmetric tensor. 


Appendix D. Low-level implementations of the SO^ li¬ 
brary 

In this section we list and describe all functions in the header 
files dlist.h, index.h, braket.h, form.h and son.h. 

dlist .h. 

• DList; ;DList(void) 

DList; :DList (const DList &L) 

DList; ;DList(int type, int i, int j) 

DList; ;DList(int type, int i) 

DList: : ~DList(void) 

The prototypes stand for the default, copy and specific 
constructors and the destructor, respectively, within the 
class DList. Note that the argument type defines the type 
of the node element. The arguments i and j are at the 
most two indices to complete the element. Pointers to the 
noList-type fields are initialised with NULL. 

• void DList : : clear () 

Deletes all nodes from DList. 

• void DList : : negate () 

Changes the sign of DList. 
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• void DList; :add(elemType) 

Adds one node at the end of DList (scans the list) and 
updates the actual pointer to be the last node. 

• void DList; :add_begin(elemType) 

Adds one node in the beginning of DList and updates the 
actual pointer to be the new hrst node. 

• void DList; :add_begin(elemType) 

Adds one node in the beginning of DList and updates the 
actual pointer to be the new hrst node. 

• void DList; ;add_end(elemType) 

Adds one node in the end of DList and updates the actual 
pointer to be the last node. 

• void DList; ; set (elemType) 

Sets data (elemType) of the node being pointed by actual 
pointer. 

• void DList ;; set_begin () 

Changes actual pointer to point at the hrst element of 
DList (beg pointer). 

• void DList ;; set_end () 

Changes actual pointer to point at the last element of DList 
(end pointer). 

• void DList ;; set_sign ( int i) 

Sets the sign of DList. 

• void DList; ;join(DList& L) 

Joins a DList to the end of the current DList (this) and 
updates the actual pointer to be the end of the hnal DList. 

• void DList; ; loop_right() 

Shifts actual pointer to next node. If actual node is the end 
node, shift to beg node. 

• void DList ;; loopjeft () 

Shifts actual pointer to previous node. If actual node is the 
beg node, shift to end node. 

• DList DList; ; rearrange 0 

Creates and returns a new DList by copying nodes in 
DList ordered by type. The nodes that hrst appear in the 
new ordered DList are d’s (type=2) and then all other ele¬ 
ments: b (type=0) and h* (type=l) unordered. Constant 
elements are removed. 

• void DList; ; remove (unsigned int type) 

Removes the hrst element with data.get_type()== 
type found in DList. Updates actual pointer to be the 
hrst node. 


• void DList ; ; remove_actual () 

Removes the element for which the actual pointer, actual, 
is pointing at in DList. 

• void DList; ; shift_right 0 

Shifts actual pointer to next node. If actual node is the end 
node, stops. 

• void DList; ; shiftjeft () 

Shifts actual pointer to previous node. If actual node if 
hrst node (begin), stops. 

• void DList ; ; swap_next () 

Swaps the actual node with the next node of DList. 

• elemType DList; ;get() 

Returns elemtype of the node being pointed by actual (cur¬ 
rent element). 

• int DList; ;getSign() 

Returns the sign of DList. 

• vector<int> DList; ;getlds() 

Creates and returns an integer vector sequence container 
with the ids (data helds) of h’s and h' ’s elements. 

• void DList; ; getBandBdagger Ids (bool feBandBd, 
vector <string> feidO, vector<string> feidl, int 
fesign) 

Updates integer vector sequence containers idO and idl 
with ids (data helds) of b’s and b^’s elements, respectively; 
the parameter sign is updated with the sign of DList and 
BandBd is a boolean which is true if DList contains at 
least one b or h’ and false otherwise. 

• void DList; ; getDeltalds (bool feAllDeltas, vector 
<string> &idO, vector<string> feidl , int fesign) 

Updates integer vector sequence containers idO and idl 
with the hrst and second ids (data helds) of 6 elements, 
respectively; the parameter sign is updated with the sign 
of DList and AllDeltas is a boolean which is true if all 
elements in DList are of 6 type and false otherwise. 

• void DList; ; getBandBdaggerAndDeltasIds ( vector 

<string> &idO, vector<string> feidl , vector< 
string> &id2, vector<string> &id3, int fesign) 

Updates integer vector sequence containers idO and idl 
with ids (data helds) of b’s and h’ ’s elements, respec¬ 
tively; updates integer vector sequence containers id2 
and ids with hrst and second data helds of d’s elements, 
respectively. The parameter sign is updated with the sign 
of DList. 

• int DList; ;numDeltas() 

Returns the number of elements of type 6 (type=2). 


• int DList; :numBs() 

Returns the number of elements of type b (type=0). 

• bool DList: ; searchjast (unsigned int typel) 

Search the last element with data.get_type()==typel 
found in DList. Returns true a node was found. 

• bool DList:: search_flrst (unsigned int typel) 

Search the first element with data.get_type()==type 
found in DList and returns true if the symbol is not the 
first element, and false otherwise. 

• bool DList:: searcb_first (unsigned int typeO, 
unsigned int typel) 

Checks for the order of appearance in DList of types 
typeO and typel. Returns true if order of appearance is 
the same as the parameter’s order and false otherwise. 

• bool DList: :search_elem( unsigned int typel) 

Searches for the element with data. get_type () ==type 
found in DList. Returns true if the element is found, false 
otherwise. 

• bool DList: : check 0 

Verifies if the number of b’s and ’s matches and whether 
the number each one is less or equal than N of SO{2N). 

• bool DList: :checkDeltaIndex() 

Checks the indexes of 6 elements and if the indexes in 6 are 
equal. Returns true if each 6 is not zero, false otherwise. 

• bool DList: :check_num() 

Verifies if the number of b’s and b^’s is less or equal than 
N of SO(2V) if so the function returns true if not returns 
false. 

• bool DList :: check_same_num () 

Verifies if the number of b’s and b^’s matchs; returns true 
if they match and false otherwise. 

• bool DList :: isActualLast () 

Returns true if actual pointer is pointing to the last (end) 
node of DList. 

• bool DList : : isEmpty () 

Returns true if DList has no nodes. 

• bool DList: :hasNoDeltas() 

Returns true if there is no elements of type 6 in DList. 

• bool DList: :hasOnlyDeltas() 

Returns true if all nodes in DList are of 6 type. 

• bool DList: :hasRepeatedIndex() 

Returns true if there is elements with the same id (data 
fields) in the DList (repeated ids). 
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• DListfe DList: :operator=(const DListfe L) 

Copies a DList. 

• const DList DList: : operator- () const 

Negates operator, change sign of DList. 

• friend DList* copy (DList *L) 

Creates and returns a pointer to a new copy of a DList. 

• friend DList contract_deltas (DList &L, bool 

braketmode) 

Applies the following identity bi *b\ - ^ij - b'j * bi, input 
DList&L keeps delta term and the function returns the 
swapped term; the parameter braketmode- if true and if 
last element in DList is a b, then the L is cleared. Returns 
expression with bi * b\ swapped or empty expression if b 
is the last term in L. 

• friend DList ordering (DList &L, bool braketmode 

) 

Order only the / 7 ’s (to the left hand side) and Z?* ’s (to the 
right hand side) terms. Applies the following identity 
b\ * bj = dij - bi * b\, input DList L keeps delta term 
and function returns the swapped term, braketmode 
- if true and if last element in DList is a b, then L is 
cleared. Returns expression with b \ * bi swapped or empty 
expression if b is the last term in L. 

• friend string printDeltas (DList &L) 

Creates and returns a string with the deltas and constants 
of a DList. 

• friend DListfe operator* (DListfe L, elemType j) 

Adds element j to the end of DList and returns a pointer 

to DList. 

• friend DList operator* (const DListfe L, const 
DList feM) 

Creates and returns a new DList that joins two DLists by 
order of parameters. 

• friend DListfe operator-(DListfe L, elemType j) 

Negates the sign of DList and adds the elemtype j at 
end of it. 

• friend DListfe operator , (DListfe L, elemType j) 

Adds element j to the end of DList and returns a pointer 

to DList. 

• friend ostreeunfe operator<< (ostreeunfe out, DList 

feL) 

Sends to output streaun ostreaun a string with the 
corresponding expression of the DList. 


• friend ostreamfe operator<< (ostreamfe out, DList 

*L) 

Sends to output stream ostream a string with the 
corresponding expression of the DList. 

• friend DListfe operator<< (DList &L, elemType j) 

Adds element j to the end of DList and returns a pointer 

to DList. 

• friend DListfe operator « (DList &L, DList &M) 

Copies DList; creates and returns a new DList with nodes 
of both old and new DLists. Sign is the product of both 
products. 

• friend bool operator== (DList &L, DList &M) 

Returns true if two DLists are equal. 

index, h. 

• int newldxdnt i) 

See dehnition in Section 3.2. 

• int newldx( string i) 

Stores a new index of type string . 

• void newld (string i) 

Stores a new index of type string . 

• string getidxdnt i) 

Returns the index placed at the position i. 

• int Idx_size() 

Returns index list size. 

• string IndexList () 

Returns index list in string of the form "Indices 

• void printIDXO 

Prints index list. 

• string makeld (string a, int id) 

Returns a+id in a string format. 

• template < class T> string ToStringd number) 

Converts to string. 

braket . h. 

• void setSimplifyIndexSum () 

See dehnition in Section 3.2. 

• void unsetSimplifyIndexSumO 

See dehnition in Section 3.2. 

23 


• BraketOneTerm ; ;BraketOneTerm() 

Constructor. 

• BraketOneTerm; ; BraketOneTerm (const DList fed) 

Constructor without constant part and index zero; d is a 
DList expression. 

• BraketOneTerm; ; BraketOneTerm (int indexin, 
string constpartin, const DList fed) 

Constructor; indexin is the index of the expression, 
constpartin is the constant part and d is a DList 
expression. 

• BraketOneTerm; ; BraketOneTerm (int indexin, 
string constpartin, list<DList> termin) 

Constructor; indexin is the index of the expression, 
constpartin is the constant part and termin is a 
list<DList> expression. 

• BraketOneTerm; ; BraketOneTerm (int indexin, 
string constpartin, BraketOneTermfe termin) 

Constructor; indexin is the index of the expression, 
constpartin is the constant part and termin is a 
BraketOneTerm expression. 

• BraketOneTerm ; ; ^BraketOneTerm () 

Destructor; clears all allocated memory. 

• void BraketOneTerm; ; clear () 

Clears all allocated memory and sets default parameters. 

• 1 i St <DList>& BraketOneTerm ; ; GetTerm () 

Returns and sets the term part. 

• string & BraketOneTerm ; ; GetConst () 

Returns (and sets) the constant part. 

• int& BraketOneTerm ;; Getindex () 

Returns (and sets) the index sum part. 

• bool BraketOneTerm; ; Simplify (QPMode oper) 
Simplihes current expression term. 

oper term mode (bra, braket, ket or none) 

return true if expression term is empty, false otherwise 

• bool BraketOneTerm; ; checkindex () 

Checks global index in expression term; returns true if |in- 
dex| is equal to 0 or A of SO(2A), otherwise returns false. 

• void BraketOneTerm; ;expfromForm( string a) 

To pass an expression from form. 


• void BraketOneTerm ; ; rearrange 0 

Orders nodes of DList in Braket. First deltas and then h’s 
and h ’ ’s, and removes the identity node when 6, b or h ’ are 
present. 

• bool BraketOneTerm; ;isempty() 

Returns true if expression is empty. 

• bool BraketOneTerm; ;EvaluateToDeltas(OPMode 

oper) 

Evaluates the expression to deltas; the mode of oper can 
be bra, braket, ket or none. This function returns true if 
the term is empty (or gives zero) otherwise returns false. 

• bool BraketOneTerm ; ; EvaluateToLeviCivita ( 

OPMode oper) 

Evaluates the expression to Levi-Civita tensors with 
eventual d’s; the mode of oper can be bra, braket, ket or 
none. This function returns true if the term is empty (or 
gives zero) otherwise returns false. 

• void BraketOneTerm; ;neg() 

Negates BraketOneTerm. 

• BraketOneTerm BraketOneTerm; ; operator*(const 
string constval) 

Overloads operator for BraketOneTerm * constval. 

• BraketOneTerm BraketOneTerm; ;operator*= (const 

string constval) 

Overloads operator for BraketOneTerm *= constval. 

• BraketOneTerm BraketOneTerm; ; operator*(const 
BraketOneTerm &L) 

Overloads operator for BraketOneTerm * L. 

• BraketOneTerm BraketOneTerm; ;operator*= (const 

BraketOneTerm &L) 

Overload operator for BraketOneTerm *= L. 

• friend BraketOneTerm operator-(const 
BraketOneTerm &L) 

Negates operator. 

• friend ostreamfe operator<< (ostreamfe out, const 

BraketOneTerm &L) 

Stream operator. 

• Braket; ; Braket (void) 

Constructor; the default expression mode is none. 

• Braket; ;Braket(const DList &d0) 

Constructor without constant part and index zero; dO is 
DList expression. 
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• Braket; ;Braket(int id, string a, DList dO) 

Constructor; the default expression mode is none and id 
is the index of the expression, a is the constant part and dO 
is the DList expression. 

• Braket; ; Braket(int id, string a, DList dO, 

OPMode op) 

Constructor, default expression mode is none. 

id - index of the expression 

a - constant part 

dO - DList expression 

op - Braket type, i.e., bra/ket/braket/none 

• Braket; ; Braket (const Braket &L) 

Constructor. 

L - Braket expression 

• Braket; ; Braket (int id, string a, const Braket &L 

, OPMode op) 

Constructor. 

id - index of the expression 
a - constant part 

L - Braket expression, ignores current constant part of L 
op - Braket type, i.e., bra/ket/braket/none 

• Braket; ; Braket (BraketOneTerm term) 

Constructor, default expression mode is none, 
term - BraketOneTerm expression 

• Braket; ; Braket (BraketOneTerm term, OPMode op 

) 

Constructor. 

term - BraketOneTerm expression 
op - Braket type, i.e., bra/ket/braket/none 

• Braket; ;~Braket() 

Destructor, clears all allocated memory. 

• void Braket; ; clear 0 

Clears all allocated memory and sets default parameters. 

• void Braket; ;expfromForm(vector<string> a) 

To pass an expression from form. 

• OPModefe Braket; ;Type() 

Returns the current expression type, it also allows to set 
new expression type. Expression types: bra/ket/braket or 
none. 

• void Braket ; ; mode () 

Prints the Braket expression mode, ie, the type of Braket; 
bra/ket/braket or none. 


• void Braket: : evaluate (bool onlydeltas=true) 

See Section 3.2. 

• void Braket : : simplify () 

Simplifies expression. Applies the following rules: 
bj |0) = 0 and (0| = 0. In (0|... |0) the number of fi,- must 

be equal to the number of b\. It also checks for numeric 
deltas and evaluates them. 

• void Braket: : rearrange 0 

Orders nodes of DList in Braket. First deltas and then h’s 
and h^’s and removes the identity node when deltas, b or 
h’ are present. 

• void Braket: :checkDeltaIndex() 

Checks index in the deltas. 

• void Braket : : gindexsetnull () 

Sets to zero the index sum of each expression term. 

• void Braket: : checkindex () 

Checks global index in expression 
term if setSimplifyIndexSum () or 

FlagSimplifyGlobalIndexSum () is active, returns 
true if |index| is equal to 0 or A of SO(2N), otherwise 
returns false. 

• void Braket: :setON() 

Activates expression term numbering for output writing 
for each term Local R?=. 

• void Braket : : setOFF () 

Deactivates expression term numbering for output writing 
for each term Local R?=. 

• int Braket :: size () 

Returns number of terms in current expression. 

• BraketOneTermfe Braket :: Get (int pos) 

Returns expression term at position given by pos. 

• intfe Braket: : Getlndex( int pos) 

Retums/sets the index sum of the term given by pos. 

• Braket Braket: : operator =( const Braket &L) 
Overloads operator for Braket = L. 

• Braket Braket: : operator +( const Braket &L) 
Overloads operator for Braket + L. 

• Braket Braket: : operator +=( const Braket &L) 
Overloads operator for Braket += L. 

• Braket Braket: : operator-(const Braket &L) 
Overloads operator for Braket - L. 


• Braket Braket: :operator-= (const Braket &L) 

Overloads operator for Braket -= L. 

• Braket Braket: : operator*(const Braket &L) 

Overloads operator for Braket * L. 

• Braketfe Braket: :operator*= (const Braket &L) 
Overloads operator for Braket *= L. 

• Braket Braket: : operator*(const string constval) 

Overloads operator for Braket * constval, i.e., the con¬ 
stant part. 

• Braket Braket: :operator*= (const string constval 

) 

Overloads operator for Braket *= constval, i.e., the 
constant part. 

• friend Braket operator-(const Braket &L) 

Overloads operator for negate, -L. 

• friend OPMode operator*(const OPMode a, const 
OPMode b) 

Calculates the mode for the multiplication. Returns mode 
of the multiplication, if this results in an invalid operation 
the program exit, a is the mode of the left operand, b 
is the mode of the right operand. 

• friend OPMode operator*(const OPMode a, const 
OPMode b) 

Calculates the mode for the sum. Returns mode of the 
sum, if this results in an invalid operation the program 
exit, a is mode of the left opercind e b is mode of the 
right operand. 

• friend OPMode operator-(const OPMode a, const 
OPMode b) 

Calculates the mode for the subtraction. Returns mode of 
the subtraction, if this results in an invalid operation the 
program exit, a is the mode of the left operctnd e b is 
mode of the right operctnd. 

• friend ostreamfe operator<< (ostreamfe out, const 

Braket &L) 

Writes expression to ostream. 

• friend string & operator «( string & out, const 
Braket &L) 

Writes expression to string . 

• friend string & operator* ( string & out, const 
Braket &L) 

Writes expression to string . 
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• ostreamfe operator<<(ostream& out, const 
OPMode &a) 

Gets the mode of the expression a of the current 
expression and returns the mode in ostreeim. 

form. h. 

• void setFormRenumberO 
void unsetFormRenumberO 
void setFormIndexSumO 
void unsetFormIndexSumO 

See Section 3.2. 

• string FormField (const string fieldname, const 

nnsigned int numUpperlds , const nnsigned int 
numLowerlds, const FuncProp funcp) 

See Section 3.2. 

• void CallFormCBraket feexp, bool print=true, 
bool all=trne, string newidlabel="j") 

See Section 3.2. 

• ToForm: : ToForm(void) 

Constructor. 

• ToForm; : ~ToForm() 

Destructor. 

• void ToForm; ; clear () 

Clears all allocated memory and sets default values. 

• bool ToForm; ; function (string f) 

Stores a field name. 

• void ToForm; ; contractions (string f) 

Stores all field contractions. 

• ToForm; ; string getFCO 

Returns all type of contractions for the fields. 

• string ToForm; ; getFunction () 

Returns all the field names. 

• void ToForm; ;setFilename (string name) 

Sets the beginning of a input/output FORM file. 

• string ToForm; ;flle() 

Returns the beginning of a input/output FORM file. 

• string & ToForm; ;rpath() 

Returns the full path name and form binary file 

• void ToForm; ;run(Braket feexp, bool print, 
bool all, string newidlabel) 

Simplify expression in FORM. Creates file input for 


FORM, runs the FORM program and returns the result to 
file and/or screen. 

exp - Braket expression to be simplified in FORM, the 
result is written back. 

print - if TRUE prints final result to screen 

all - if TRUE writes all the expression members sepa¬ 
rately in ouput EORM file, if EALSE only writes the full 
result together 

newidlabel - label to be used when the option to sum 
index is active 

• bool ToForm; ; getIndexSum () 

Returns the state of the indexSum flag. 

• void ToForm; ;setIndexSum (bool flag) 

Sets the state of the indexSum flag . 

• void ToForm; ;setRenumber(bool flag=true) 

Sets "renumber 1;" in FORM input file. This option is used 
to renumber index in order to allow further simplifications. 
However, in big expressions this must be avoid since it in¬ 
creases the computational time in FORM. The best way to 
use is simplify the expression with FORM with this op¬ 
tion unset, and then send a second time to FORM with this 
option active. By default this option is unset. 

• bool ToForm; ; getRenumberOption () 

Returns the state of the f ormRenumber flag. 

• ToFormfe ToForm; ; operator«(const string fefunc) 
Overloads operator for ToForm << func. 

• ToFormfe ToForm; ;operator+(const string fefunc) 
Overloads operator for ToForm + func. 

son. h. 

• void setDimdnt n) 
int getDimO 

void CleanGlobalDecK) 

void setVerbosity ( Verbosity verb) 

Verbosity getVerbosity () 

Braket Bop(std; ; string startid="i") 

Braket BopIdnumO 

The description of these functions was done in Section 3.2. 

• bool GroupEvenO 

Returns Group parity. 

• ostreaunfe operator« (ostreaunfe out, const 
Verbosity fea) 

Returns current ostream verbosity level. Returns ostream 
for output. 
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• size_t getCurrentRSS () 

Returns the current resident set size (physical memory use) 
measured in bytes, or zero if the value cannot be deter¬ 
mined on this OS. 

• size_t getPeakRSSO 

Returns the peak (maximum so far) resident set size (phys¬ 
ical memory use) measured in bytes, or zero if the value 
cannot be determined on this OS. 

• void print_process_mem_usage() 

Prints memory stats (current and peak resident set size) in 
MB. 

In addition, we have also implemented C-i-H- macros that sim¬ 
plifies the call for the functions; its description was done in Sec¬ 
tion 3.2. Notice that the stringising macro operator #a causes 
the argument to be enclosed in double quotation marks. 

• bb(id) 

This makes the call DListCO, newldx(id)), id is the 
index in string format or enclosed in quotation marks. 

• bbt(id) 

This makes the call DListCl, newldx(id) ), id is the 
index in string format or enclosed in quotation marks. 

• b(id) 

This makes the call DListCO, newldx(#id) ), id index 
does not need to be enclosed in quotation marks. 

• bt(id) 

This makes the call DListCl , newIdxC#id) ), id index 
does not need to be enclosed in quotation marks. 

• delta Cidl, id2) 

This makes the call DListC2, newIdxC#idl) , newldx 
C#id2)) ; idl and id2 indices do not need to be enclosed 
in quotation marks. 

• identity 

This macro is a shortcut for the object DListCS , 0). 

• braCi,a,b) 

This makes the call BraketCi, #a, b, bra), a index 
does not need to be enclosed in quotation marks. 

• ketCi,a,b) 

This makes the call BraketCi, #a, b, ket), i index 
sum, a constant part without quotation marks, b DList ex¬ 
pression. 

• braketCi,a,b) 

This makes the call BraketCi, #a, b, braket), i in¬ 
dex sum, a constant part without quotation marks, b DList 
expression 


• freeCi,a,b) 

This makes the call BraketCi , #a, b, none), i index 
sum, a constant part without quotation marks, b DList ex¬ 
pression. 

• Field Ca, b, c, d) 

This makes the call ForniFieldC#a, b, c, d), a held 
name without quotation marks, b number of upper index, 
c number of lower index, d held with/without havor in¬ 
dex and symmetric/asymmetric held, returns held name in 
string format. 
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